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METHODS AND APPARATUS FOR 
DNA SEQUENCING AND DNA IDENTinCATION 

Field of the Invendon 

This invention relates in general to methods and apparatus for 
nucleic acid analysis, and, in particular to, methods and ayiparatus for DNA 
sequencing. 

Background 

The rate of determining the sequence of the four nucleotides in 
DNA samples is a major technical obstacle for further advancement of 
molecular biology, medicine, and biotechnology. Nucleic acid sequencing 
methods which involve separation of DNA molecules in a gel have beerrin use 
since 1978. The only other proven method for sequencing nucleic acids is 
sequencing by hybridization (SBH), 

The array-based approach of SBH does not require single base 
resolution in separation, degradation, synthesis or imaging of a DNA 
molecule. In the most commonly discussed variation of this method, using 
mismatch discriminative hybridization of short oligonucleotides K bases in 
length, lists of constituent K-mer oligonucleotides may be determined for 
target DNA. The sequence may be assembled through uniquely overlapping 
scored oligonucleotides. 

In SBH sequence assembly, K -1 oligonucleotides which occur 
rq)eatedly in analyzed DNA fragments due to chance or biological reasons 
may be subject to special consideration. If there is no additional information, 
relatively small fragments of DNA may be fully assembled in as much as 
every base pair (bp) is read several times. In assembly of relatively longer 
fragments, ambiguities may arise due to repeated occurrence of a K -1 
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nucleodde. This problem does not exist if mutated or similar sequences have 
to be determined. Knowledge of one sequence may be used as a template to 
correctly assemble a similar one. 

There are several approaches for sequencing by hybridiiation. 
5 In SBH Format 1, DNA samples are arrayed and labelled probes are 
hybridized with the samples. Replica membranes with the same sets of 
sample DNAs may be used for parallel scoring of several probes and/or 
probes may be multiplexed. Arraying and hybridization of DNA samples on 
the nylon membranes are well developed. Each array may be reused many 
10 times. Format 1 is especially efficient for batch processing large numbers of 
sanqiles. 

In SBH Format 2, probes are arrayed and a labelled DNA 
sample fragment is hybridized to the arrayed probes. In this case, the 
complete sequence of one fragment may be determined from simultaneous 

15 hybridization reactions with the arrayed probes. For sequencing other DNA 
fragments, the same oligonucleotide array may be reused. The arrays may be 
produced by spotting or in situ variant of Format 2, DNA anchors are arrayed 
and ligation is used to determine oligosequences present synthesis. Specific 
hybridization has been demonstrated. In a variant of Format 2, DNA anchors 

20 are arrayed and ligation is used to determine oligosequences present at the end 
of target DNA. 

In Format 3, two sets of probes are used. One set may be in 
the form of arrays and another, labelled set is stored in multiwell plates. In 
this case, target DNA need not be labelled. Target DNA and one labelled 
25 probe are added to the arrayed set of probes. If one attached probe and one 
labelled probe both hybridize contiguously on the target DNA, they are 
covalently ligated, producing a sequence twice as long to be scored. The 
process allows for sequencing long DNA fragments, e.g. a complete bacterial 
genome, without DNA subcloning in smaller pieces. 
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In the present invention, SBH is ^lied to the efficient 
identification and sequencing one or more DNA samples in a short period of 
time. The procedure has many applications in DNA diagnostics, forensics, 
and gene mapping. It also may be used to idoitify mutations responsible for 
S genetic disorders and other traits, to assess biodiversity and to produce many 
other types of data dq)mdem on DNA sequence. 

Summary of the Invention 

As mentioned above, Formal I SBH is appropriate for the 
simultaneous analysis of a large set of samples. Parallel scoring of thousands 
10 of samples on large arrays may be applied to one or a few samples are in 
thousands of independent hybridization reactions using small pieces of 
membranes. The identification of DNA may involve 1-20 probes and the 
identification of mutations may in some cases involve more than 1000 probes 
specifically selected or designed for each sample. For identification of the 
15 nature of the mutated DNA segments, specific probes may be synthesized or 
selected for each mutation detected in the first round of hybridizations. 

According to the present invention, DNA samples may be 
prepared in small arrays which may be separated by appropriate spacers, and 
which may be simultaneously tested with probes selected from a set of 
20 oligonucleotides kept in multiwell plates. Small arrays may consist of one or 
more samples. DNA samples in each small array may consist of mutants or 
individual samples of a sequence. Consecutive small arrays which form larger 
arrays may represent either replication of the same array or samples of a 
di£ferent DNA fiagment. A universal set of probes consists of sufficient 
25 probes to analyze any DNA fragment with prespecified precision, e.g. with 
Tcspea to the redundancy of reading each bp. These sets may include more 
probes than are necessary for one specific fragment, but fewer than are 
necessary for testing thousands of DNA samples of different sequence. 
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DNA or allele identification and a diagnostic sequencing process 
may include the steps of: 



1) Selection of a subset of probes from a dedicated, lepresentative or 
universal set to be hybridized with each of a plurality small arrays; 
S 2) Adding a first probe to each subarray on each of the arrays to be 
analyzed in parallel; 

3) Performing hybridization and scoring of the hybridization results; 

4) Stripping off previously used probes and repeating remaining probes 
that are to be scored; 

10 5) Processing the obtained results to obtain a final analysis or to 
determine additional probes to be hybridized; 

6) Performing additional hybridizations for certain subarrays; and 

7) Processing complete sets of data and computing obtaining a final 
analyns. 



15 The present invention solves problems in fast identification and 

sequencing of a small number of nucleic acid samples of one type (e.g. DNA, 
RNA) and in parallel analysis of many sample types by using a presynthesized 
set of probes of manageable size and samples attached to a support in the form 
of subarrays. Two approaches have been combined to produce an efficient 

20 and versatile process for the determination of DNA identity, for DNA 
diagnostics, and for identification of mutations. For the identification of 
known sequences a small set of shorter probes may be used in place of a 
longer unique probe. In this case, there may be more probes to be scored, but 
a universal set of probes may be synthesized to cover any type of sequence. 

25 For example, a full set of 6-mers or 7-mers are only 4,096 and 16,384 
probes, respectively. 

Full sequencing of a DNA fragment may involve two levels. 
One level is hybridization of a sufficient set of probes that cover every base 
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at least once. For this purpose^ a specific set of probes may be synthesized 
for a standard sample. This hybridization data reveals whether and where 
mutations (differences) occur in non-standard samples. To determine the 
identity of the changes, additional specific probes may be hybridized to the 
5 sample. In anoth^ embodiment, all probes from a universal set may be 
scored. 

A universal set of probes allows scoring of a relatively small 
number of probes per sample in a two step process without unacceptable 
expenditure of time. The hybridization process involves successive probings, 
10 in a first step of computing an optimal subset of probes to be hybridized first 
and, then, on the basis of the obtained results, a second step of determining 
additional probes to be scored from among those in the existing universal set. 

The use of an array of sample arrays avoids consecutive scoring 
of many oligonucleotides on a single sample or on a small set of samples. 
15 This approach allows the scoring of more probes in parallel by manipulation 
of only one physical object. By combining the use of the subarray formed 
with the universal set of probes and the four step hybridization process, a 
DNA sample 1000 bp in length may be sequenced in a relatively short period 
of time. If the sample is spotted at 50 subarrays in an array and the array is 
20 reprobed 10 times, 500 probes may be scored. This number of probes is 
highly sufficient. In screening for the occurrence of a mutation, 
approximately 335 probes may be used to cover each base three times. If a 
mutation is present, several covering probes will be affected. These negative 
probes may map the mutation with a two base precision. To solve a single 
25 base mutation noapped with this precision, an additional 15 probes may be 
employed. These probes cover any base combination for the two questionable 
positions (assuming that deletions and insertions are not involved). These 
probes may be scored in one cycle on SO sidiairays vAAch contain the given 
sample. In the implementation of a multiple label color scheme 
30 (multiplexing), two to six probes labelled with different fluorescent dyes may 
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be used as a pool, thereby reducing the number of hybridization cycles and 
shortoiing the sequencing process. 

In more complicated cases, there may be two close mutations 
or insertions. They may be handled with more probes. For example, a ttiree 
5 base insertion may be solved with 64 probes. The most complicated cases 
may be approached by several steps of hybridization, and the selecting of a 
new set of probes on the basis of results of previous hybridizations. 

If subanays consists of tens or hundreds of samples of one 
type, then several of them may be found to contain one or more changes 

10 (mutations, insertions, or deletions). For each segment where mutation 
occurs, a specific set of probes may be scored. The total number of probes to 
be scored for a type of sample may be several hundreds. The scoring of 
replica arrays in parallel allow scoring of hundreds of probes in a relatively 
small number of cycles. In addition, compatible probes may be pooled. 

15 Positive hybridizations may be assigned to the probes selected to check 
particular DNA segments because these segments usually differ in 75% of 
their constituent bases. 

By using a larger set of longer probes, longer targets may be 
conveniently analyzed. These targets may represent pools of shorter 

20 fragments sudi as pools of exon clones. 

The multiple step approach, which minimizes the number of 
necessary probes, may employ a specific hybridization scoring method to 
define the presence of heterozygotes (sequence variants) in a genomic segment 
to be sequenced from a diploid chromosomal set. There are two possibilities: 

25 i) the sequence from one chromosome represents a basic type and the 
sequence from the other represents a new variant; or, ii) both chromosomes 
contain new, but different variants. In the first case, the scanning step 
designed to map changes gives a maximal signal difference of two-fold at the 
heterozygotic position. In the second case, there is no masking; only a more 
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complicated selection of the probes for the subsequent rounds of hybridizations 
may be required. 

Scoring two-fold signal differences required in the first case 
may be achieved efficiently by comparing conesponding signals with controls 
5 containing only the basic sequence type and with the signals from other 
analyzed samples. This approach allows determination of a relative reduction 
in the hybridization signal for each particular probe in the given sample. This 
is significant because hybridization efficiency may vary more than two-fold 
for a particular probe hybridized with different DN A fragments having its fiill 
10 match target. In addition, heterozygotic sites may affect more than one probe 
depending on the number of oligonucleotide probes. Decrease of the signal 
for two to four consecutive probes produces a more significant indication of 
heterozygotic sites. The leads may be checked by small sets of selected 
probes among which one or few probes are suppose to give full match signal 
IS which is on average eight-fold stronger than the signals coming from 
mismatch-containing duplexes. 

Partitioned membranes allow a very flexible organization of 
experiments to accommodate relatively larger numbers of samples representing 
a given sequence type, or many different types of samples represented with 
20 smaller number of samples. A range of 4>2S6 samples can .be handled with 
particular efficiency. Subarrays within this range of numbers of dots may be 
designed to match the configuration and size of standard multiwell plates used 
for storing and labelling oligonucleotides. The size of the subarrays may be 
adjusted for different number of samples, or a few standard subarray sizes 
25 may be used. If all samples of one ^pe do not fit in one subarray, additional 
subarrays or membranes may be used and processed with the same probes. 
In addition, by adjusting the number of replicas for each subarray, the dme 
for completion of tdentificadon or sequencing process may be varied. 
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DetaUed Description 

Example I 
Prepaxadon of a Universal Set of Probes 

Two types of universal sets of probes may be prepared. The 
5 first is a ccmiplete set (or at least a noncomplementary subset) of relatively 
short probes. For example, all 4096 (or about 2CXX) non-complementary) 

6- mm, or all 16,384 (or about 8,000 non-complementary) 7-mers. Full 
noncomplementary subsets of S-mers and.longer probes are less convenient in 
as much as they include 32,000 or more probes. 

10 A second type of probe set is selected as a small subset of 

probes still sufficient for reading every bp in any sequence with at least with 
one probe. For example, 12 of 16 dimers are sufficient. A small subset for 

7- mers/8-mer and 9-mers for sequencing double stranded DNA may be about 
3000, 10,000 and 30,000 probes, respectively. 

15 Probes may be prepared using standard chemistry with one to 

three non-specified (mixed A,T,C and G) or universal (e.g. M base» inosine) 
bases at the ends. If radiolabelling is used, probes may have an OH group at 
the 5' end for kinasing by radiolabeled phosphorous groups. Alternatively, 
probes labelled with fluorescent dyes may be employed. Other types of 

20 probes like PNA (Protein Nucleic Acids)or probes containing modified bases 
which change duplex stability also may be used. 

Probes may be stored in barcoded multiwell plates. For small 
numbers of probes, 96-weli plates may be used; for 10,000 or more probes, 
storage in 384* or 864-well plates is preferred. Stacks of S to SO plates are 

25 enough to store all probes. Approximately S pg of a probe may be sufficient 
for hybridization with one DNA sample. Thus, from a small synthesis of 
about SO §1$ per probe, ten million samples may be analyzed. If each probe 
is used for every third sample, and if each sample is KXX) bp in length, then 

SDBSnnnESSRn(BllL£2$) 
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over 30 billion bases (10 human genomes) may be sequenced by a set of S»000 
probes. 

Example 2 
Preparation of DNA Samples 

5 DNA fragments may be prq>ared as clones in M13, plasmid or 

lambda vectors and/or prepared directly from genomic DNA or cDNA by 
PCR or other amplification methods. Samples may be prepared or di^)ehsed 
in multiwell plates. About 100-1000 ng of DNA samples may be prqwed in 
2*500 fd of final volume. 

10 Example 3 

Preparation of DNA Arrays 

Arrays may be prepared by spotting DNA samples on a support 
such as a nylon membrane. Spotting may be performed by using arrays of 
metal pins (the positions of which correspond to an array of wells in a 

15 microliter plate) to repeated by transfer of about 20 nl of a DNA solution to 
a nylon membrane. By offset printing, a density of dots higher than the 
density of the wells is achieved. One to 25 dots may be accommodated in 1 
mih^ depending on the type of label used. By avoiding spotting in some 
preselected number of rows and columns, separate subsets (subarrays) may be 

20 formed. Samples in one subarray may be the same genomic segment of DNA 
(or the same gene) from different individuals, or may be different, overlapped 
genomic clones. Each of the subarrays may represent replica spotting of the 
same samples. In one example, one gene segment may be amplified ftom 64 
patients. For each patient, the ampliHed gene segment may be in one 96-well 

25 plate (all 96 wells containing the same sample). A plate for each of the 64 
patients is prepared. By using a 96-pin device all samples may be spotted on 
one 8 X 12 cm membrane. Subarrays may contain 64 samples* one from each 
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padent. Where the 96 subanays are identical, the dot span may be 1 mm^ and 
there may be a 1 mm space between subanays. 

Another aiq>rDach is to use membranes or plates (available from 
NUNC, N2q)erville, Illinois) which may be partitioned by physical spacers e.g. 
a plastic grid nu>lded over the membrane, the grid being similar to the sort of 
membrane SQyplied to the bottom of multiwell plates, or hydrophobic strips. 
A fixed physical spacer is not preferred for imaging by exposure to flat 
phosphor-storage screens or x-ray films. 

Example 4 
Selection and Labelling of Probes 

When an array of subanays is produced, the sets of probes to 
be hybridized in each of the hybridization cycles on each of the subarrays is 
defined. For the samples in- Example 3, a set of 384 probes may be selected 
from the universal set, and 96 probings may be performed in each of 4 cycles. 
Probes selected to be hybridized in one cycle preferably have similar G+C 
contents. 

Selected probes for each cycle are transferred to a 96-well plate 
and then are labelled by idnasing or by other labeliing procedures if they are 
not labelled (e.g. with stable fluorescent dyes) before they are stored. 

On the basis of the first round of hybridizations, a new set of 
probes may be defined for each of the subarrays for additional cycles. Some 
of the arrays may not be used in some of the cycles. For example, if only 8 
of 64 patient samples exhibit a mutation and 8 probes are scored first for each 
mutation, then all 64 probes may be scored in one cycle and 32 subarrays are 
not used. These subarrays may then be treated with hybridization buffer to 
prevent drying of the filters. 

Probes may be retrieved from the storing plates by any 
convenient apiH-oach, such as a single channel pipetting device or a robotic 
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station such as a Beckman Biomek 1000 (Beckman Instruments, FuUerton, 
California) or a Mega Two robot (Megamation, Lawrenceville, New Jersey). 
A rotwtic station may be integrated with data analysis programs and probe 
managing programs. Outputs of these programs may be inputs for one or 
more robotic stations. 

Probes may be retrieved one by one and added to subarrays 
covered by hybridization buffer. It is preferred that retrieved probes be 
placed in a new plate and labelled or mixed with hybridization buffm*. The 
preferred method of retrieval is by accessing stored plates one by one and 
pipetting (or transferring by metal pins) a sufficient amount of each selected 
probe from each plate to specific wdls in an intermediary plate. An array of 
individually addressable pipettes or pins may be used to speed up the retrieval 
process. 

Example S 
Hybridization and Scoring Process 

Labelled probes may be mixed with hybridization buffer and 
pipetted preferentially by multichannel pipettes to the subarrays. To prevent 
mixing of the probes between subarrays (if there are no hydrophilic strips or 
physical barriers imprinted in the membrane), a corresponding plastic, metal 
or ceramic grid may be firmly pressed to the membrane. Also, the volume 
of the buffer may be reduced to about 1 pti or less per mm-. The 
concentration of the probes and hybridization conditions used may be as 
described previously except that the washing buffer may be quickly poured 
over the array of subarrays to allow fast dilution of probes and thus prevent 
significant cross-hybridization. For the same reason, a minimal concentration 
of the probes may be used and hybridization time extended to the maximal 
practical level. For DNA detection and sequencing, knowledge of a "normal" 
sequence allows the use of the continuous stacking interaction phraomenon to 
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increase the signal. In addition to the labelled probe, additional unlabeUed 
probes which hybridize back to back with a labelled one may be added in the 
hybridization reaction. The amount of the hybrid may be increased several 
times. The probes may be connected by ligation. This approach may be 
5 important for resolving DNA regicms forming ^^compressions". 

In the case of radiolabelled probes, images of the filters may 
be obtained preferentially by phosphorstorage tedmology. Fluorescent labels 
may be scored by CCD cameras, confocal microscopy or otherwise. Raw 
signals are normalized based on the amount of target in each dot to properly 

10 scale and integrate data different hybridization experiments. Differences 
in the amount of target DNA per dot may be corrected for by dividing signals 
of each probe by an average signal for all probes scored on one dot. Also, 
the normalized signals may be scaled, usually from 1-100, to compare data 
from different experiments. Also, in each subanay, several control DNAs 

15 may be used to determine an average background signal in those samples 
which do not contain a full match target. Furthermore, for samples obtained 
from diploid (polyploid) scores, homozygotic controls may be used to allow 
recognition of heterozygotes in the samples. 

Example 6 

20 Diagnostics - Scoring Known Mutations or Full Gene Resequencing 

A simple case is to discover whether some known mutations 
occur in a DNA segment. Less than 12 probes may suffice for this purpose, 
for example, S probes positive for one allele, S positive for the other, and 2 
negative for both. Because of the small number of probes to be scored per 
25 sample, large numbers of samples may be analyzed in parallel. For example, 
with 12 probes in 3 hybridization cycles, 96 different genomic loci or gene 
segments from 64 patient may be analyzed on one 6 x 9 in membrane 
containing 12 x 24 subanays each with 64 dots representing the same DNA 
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segment from 64 patients. In this ^cample, samples may be prepared in 
sixty-four 96-weH plates. Each plate may represent one patient, and each well 
may represent one of the DNA segments to be analyzed. The samples from 
64 plates may be spotted in four rq>licas as four quarters of the same 
5 membrane. 

A set of 12 probes may be selected by single channel pipetting 
or a single pin transferring device (or by an array of mdividually controlled 
pipets or pins) for each of the 96 segments and rearranged in twelve 96-weU 
plates. Probes may be labelled if they are not prelabelled before storing, and 
10 then probes from four plates may be mixed with hybridization bu^er and 
added to the subaiiays preferentially by a 96-channel pipeting device. After 
one hybridization cycle it is possible to strip off previously used probes by 
incubating the membrane at 37"* to SSX in the preferably undiluted 
hybridization or washing buffer. 
IS The likelihood that probes positive for one allele are positive 

and probes positive for the other allele are negative may be used to determine 
which of the two allels is present. In this redundant scoring scheme, some 
level (about 10%) of errors in hybridization of each probe may be tolerated. 

An incomplete set of probes may be used for scoring most of 
20 the alleles, especially if the smaller redundancy is sufficient, e.g. one or two 
probes which prove the presence or absence in a sample of one of the two 
alleles. For example, with a set of four thousand S-mers there is a 91% 
chance of finding at least one positive probe for one of the two alleles for a 
randomly selected locus. The incomplete set of probes may be optimized to 
25 reflect G+C content and other biases in the analyzed samples. 

For full gene sequencing, genes may be amplified in an 
appropriate number of segments. For each segment, a set of probes (about 
one probe per 2-4 bases) may be selected and hybridized. These probes may 
identify whether there is a mutation anywhere in the analyzed segments. 
30 Segments (i.e., subarrays which contain these segments) where one or more 
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mutated sites are detected may be hybridized with additional probes to find the 
exact sequence at the mutated sites. If a DNA sample is tested by every 
second 6-mer» and a mutation is localized at the position that is surrounded by 
positively hybridized probes TGCAAA and TATTCC and covmd by three 
5 n^ative probes: CAAAAC, AAACTA and ACT ATT, the mutated nucleotides 
must be A and/or C occurring in the normal sequence at that position. They 
may be changed by a single base mutation, or by a one or two nucleotide 
deletion and/or insertion between bases AA» AC or CT. 

One approach is to select a probe that extends the positively 

10 hybridized probe TGCAAA for one nucleotide to the right, and which extends 
the probe TATTCC one nucleotide to the left. With these 8 probes 
(GCAAAA, GCAAAT, GCAAAC, GCAAAG and ATATTC, TTATTC, 
CTATTC, GTATTC) two questionable nucleotides are determined. 

The most likely hypothesis about the mutation may be 

15 determined. For example, A is found to be mutated to G. There are two 
solutions satisfied by these results. Either replacement of A with G is the 
only change or there is in addition to that change an insertion of some number 
of bases between newly determined G and the following C. If the result with 
bridging probes is negative these options may then be checked first by at least 

20 one bridging probe comprising the mutated position (AAGCTA) and with an 
additional 8 probes: CAAAGA, CAAAGT, CAAAGC, CAAAGG and 
ACTATT, TCTATT, CCTATT, GCTATT, I There are many other ways to 
select mutation-solving probes. 

In the case of diploid, particular comparisons of scores for the 

23 test samples and homozygotic control may be performed to identify 
heterozygotes (see above). A few consecutive probes are expected to have 
roughly twice smaller signals if the segment covered by these probes is 
mutated on one of the two chromosomes. 
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Example 7 

Identification of Genes (Mutations) Responsible for Genetic 
Disorders and Other Traits 



The sequendng process disclosed herein has a very low cost per 
5 bp. Also, using larger universal sets of longer probes (8-mers or 9-nieis), 
DNA fragments as long as 5-20 kb may be sequenced without subdoning. 
Furthermore, the speed of resequencing may be about 10 million 
bp/day/hybridization instrument. This performance allows for resequencing 
a large fraction of human genes or the human genome rq)eatedly from 
10 sdendficaily or medically interesting individuals. To resequence 50% of the 
human genes, about 100 million bp is checked. That may be done in a 
relatively shon period of time at an affordable cost. 

This enormous resequencing capability may be used in several 
ways to identify mutations and/or genes that encode for disorders or any other 
15 traits. Basically, mRNAs (which may be converted into cDNAs) from 
particular tissues or genomic DNA of patients with particular disorders may 
be used as starting materials. From both sources of DNA, separate genes or 
genomic fragments of appropriate length may be prepared either by cloning 
procedures or by in vitro amplification procedures (for example by PCR). If 
20 cloning is used, the minimal set of clones to be analyzed may be selected ftom 
the libraries before sequencing. That may be done efficiently by hybridization 
of a snjall number of probes, especially if a small number of clones longer 
than 5 kb is to be sorted. Cloning may increase the amount of hybridization 
data about two times, but does not require tens of thousands of PCR primers. 
25 In one variant of the procedure, gene or genomic fragments 

may be prepared by restriction cutting with enzymes like Hga I which cuts 
DNA in following way: GACGC(N5')/CTGCG(N10'). Protruding ends of five 
bases are different for different fragments. One enzyme produces sq^propriate 
fragments for a certain number of genes. By cutting cDNA or genomic DNA 
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with several enzymes in separate reactions, every gene of intmst may be 
excised appnq)riately. In one approach, the cut DNA is fiactionated by size. 
DNA fragments prepared in this way (and optionally treated with Exonuclease 
ni which individually removes nucleotides from the 3' end and increases 
S length and qiecificity of the ends) may be dispensed in the tubes or in 
multiwell plates. From a relatively sinall set of DNA adapters with a common 
portion and a variable protruding end of appropriate length, a pair of adapters 
may be selected for every gene fragment that needs to be amplified. These 
adapters are ligated and then PGR is performed by universal primers. From 

10 1000 adapters, a million pairs may be generated, thus a million different 
fragmoits may be q)ecifically amplified in the identical conditions with a 
universal pair of primers complementary to the common end of the adapters. 

If a DNA difference is found to be repeated in several patients, 
and that sequence change is nonsense or can change function of the 

15 corresponding protein, then the mutated gene may be responsible for the 
disorder. By analyzing a significant number of individuals with particular 
trails, functional allelic variations of particular genes could be associated by 
specific traits. 

This approach may be used to eliminate the need for very 
20 expensive genetic mapping on extensive pedigrees and has special value when 
there is no such genetic data or material. 

Example 8 

Scoring Single Nucleotide Polymorphisms in Genetic Mapping 

Techniques disclosed in this application are appropriate for an 
25 efficient identification of genomic fragments with single nucleotide 
polymorphisms (SNUPs). In 10 individuals by applying the described 
sequencing process on a large number of genomic fragments of known 
sequence that may be amplified by cloning or by in viu-o amplification, a 
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sufficient number of DNA segmrats with SNUPs may be idratified. The 
polymcnphic fragments are further used as SNUP markers. These markers are 
dther mapped previously (fen* example they represent mapped STSs) or they 
may be mapped through the screoiing procedure described bdow. 

SNUPs may be scored in every individual from relevant 
families or populations by amplifying markers and arraying them in the form 
of theanay of subarrays. Subarrays contain the same marker amplified from 
the analyzed individuals. For each marker, as in die diagnostics of known 
mutations, a set of 6 or less probes positive for one allde and 6 or less probes 
positive for the other allele may be selected and scored. From the significant 
association of one or a group of the markers with the disorder, chromosomal 
position of the responsible gene(s) may be determined. Because of the high 
throughput and low cost, thousands of markers may be scored for thousands 
of individuals. 

This amount of data allows localization of a gene at a resolution 
level of less than one million bp as well as localization of genes involved in 
polygenic diseases. Localized genes may be identified by sequencing 
particular regions from relevant normal and affected individuals to score a 
mutation(s). 

PCR is preferred for amplification of markers from genomic 
DNA. Each of the markers require a specific pair of primers. The existing 
markers may be convertible or new markers may be defined which may be 
prepared by cutting genomic DNA by Hga I type restriction enzymes, and by 
ligation with a pair of adapters as described in Example 7. 

SNUP markers can be amplified or spotted as pools to reduce 
the number of indepotdent amplification reactions. In this case, more probes 
are scored per one sample. When 4 markers are pooled and spotted on 12 
replica membranes, then 48 probes (12 per marker) may be scored in 4 
cycles. 
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ExampleP 

Detecdon and Verification of Identity of DNA Fragments 

DNA fragments generated by restriction cutting, cloning or in 
vitro amplification (e.g. PCR) frequently may be identified in a experiment. 
S Identification may be performed by verifying the presence of a DNA band of 
specific size on gd dectrophoresis. Alternatively, a specific oligonucleotide 
may be prepared and used to verify a DNA sample in question by 
hybridization. The procedure developed here allows for more efficioit 
identification of a large number of samples without preparing a specific 
10 oligonucleotide for each fragment. A set of positive and negative probes may 
be selected from the universal set for each fragment on the basis of the known 
sequences. Prol)es that are selected to be positive usually are able to form one 
or a few overlapping groups and negative probes are spread over the whole 
insert. 

15 This technology may be used for identification of STSs in the 

process of their mapping on the YAC clones. Each of the STSs may be tested 
on about 100 YAC clones or pools of YAC clones. DNAs from these 100 
reactions possibly are spotted in one subarray. Different STSs may represent 
consecutive subarrays. In several hybridization cycles, a signature may be 

20 generated for each of the DNA samples, which signature proves or disproves 
existence of the particular STS in the given YAC clone with necessary 
confidence. 

To reduce the number of independent PCR reactions or the 
number of independent samples for spotting, several STSs may be amplified 
25 simultaneously in a reaction or PCR samples may be mixed, respectively. In 
this case more probes have to be scored per one dot. The pooling of STSs is 
indqiendent of pooling YACs and may be used on single YACs or pools of 
YACs. This scheme is especially attractive when several probes labelled with 
different colors are hybridized together. 
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In addition to confirmation of the existence of a DNA iragment 
in a sample, the amount of DNA may be estimated using intensities of the 
hybridization of several separate probes or one or more pools of probes. By 
comparing obtained intensittes with intensities for control samples having a 
5 known amount of DNA, the quantity of DNA in all spotted samples is 
determined simultaneously. Because only a few probes are necessaiy for 
identification of a DNA fiagment, and thm are N possible probes that may 
be used for DNA N bases long, this application does not require a laige set 
of probes to be sufficient for identification of any DNA segment. From one 
10 diousand 8-mers, on average about 30 full matching probes may be selected 
for a 1000 bp fragment. 

Example 10 

Identification of Infectious Disease Organisms and Their Variants 

DNA-based tests for the detection of viral, bacterial, fungal and 
15 other parasitic organisms in patients are usually more reliable and less 
expensive than alternatives. The major advantage of DNA tests is to be able 
to identify specific strains and mutants, and eventually be able to apply more 
effective treatment. Two applications are described below. 

The presence of 12 known antibiotic resistance genes in 

20 bacterial infections may be tested by amplifying tiiese genes. The amplified 
products from 128 patients may be spotted in two subarrays and 24 subanays 
for 12 genes may then be repeated four times on a 8 x 12 cm membrane. For 
each gene, 12 probes may be selected for positive and negative scoring. 
Hybridizations may be performed in 3 cycles. For these tests, as for the tests 

25 in Example 9, a much smaller set of probes is most likely to be universal. 
For example, from a set of one thousand 8-mers, on average 30 probes are 
positive in 1000 bp fragments, and 10 positive probes are usually sufficient 
for a highly reliable identification. As described in Example 9, several genes 
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may be amplified and/or spotted together and the amount of the given DNA 
may be determined. The amount of amplified gene may be used as an 
indicator of the level of infecdon. 

Another example involves possible sequencing of one gene or 
S the whole genome ofan HIV virus. Because ofispiddiversification, the virus 
poses many difficulties for sdection of an optimal therapy. DNA fragments 
may be amplified from isolated viruses from up to 64 patients and resequenced 
by the described procedure. On the basis of the obtained sequence the optimal 
therapy may be selected. If there is a mixntie of two virus types of which one 

10 has the basic sequence (similar to the case of heterozygotes), the mutant may 
be identified by quantitative comparisons of its hybridization scores with 
scores of other samples, especially control samples containing the basic virus 
type only. Scores twice as small may be obtained for three to four probes that 
cover the site mutated in one of the two virus types present in the sample (see 

IS above). 



Example II 

Forensic and Parental Identification Applications 

Sequence polymorphisms make an individual genomic DNA 
unique. This permits analysis of blood or other body fluids or tissues from 

20 a crime scene and comparison with samples from criminal suspects. A 
sufficient number of polymorphic sites are scored to produce a unique 
signature of a sample. SBH may easily score single nucleotide polymorphisms 
to produce such signatures. 

A set of DNA fragments (10-1000) may be amplified from 

25 samples and suspects. DNAs from samples and suspects representing one 
fragment are spotted in one or several subarrays and each subarray may be 
replicated 4 times. In three cycles, 12 probes may determine the presence of 
allele A or B in each of the samples, including suspects, for each DNA locus. 
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Matching the pattmis of samples and suspects may lead to discovery of the 
suspect reqxmsible for the chme. 

The same procedure may be a^licable to prove or disprove Uie 
identity of parents of a child. DNA may be prepared and polymorphic loci 
amplified from the child and adults; patterns of A or B alleles may be 
determined by hybridization for each. Comparisons of the obtained patterns, 
along with positive and negative controls, aide in the determination of familial 
relationships. In this case» only a significant portion of the alleles need match 
with one parent for identification. Large numbers of scored loci allow for the 
avoidance of statistical errors in the procedure or of masking effects of de 
novo mutations. 

Example 12 

Assessing Genetic Diversity of Populations or Species and 
Biological Diversity of Ecological Niches 

Measuring the frequency of allelic variations on a significant 
number of loci (for example, several genes or entire mitochondrial DNA) 
permits development of different types of conclusions, such as conclusions 
regarding the impact of the environment on the genotypes, history and 
evolution of a population or its susceptibility to diseases or extinction, and 
others. These assessments may be performed by testing specific known alleles 
or by full resequendng of some loci to be able to define de novo mutations 
which may reveal fine variations or presence of mutagens in the environment. 

Additionally, biodiversity in the microbial world may be 
surveyed by resequendng evolutionarily conserved DNA sequences, such as 
the genes for ribosomal RNAs or genes for highly conservative proteins. 
DNA may be prepared from the environment and particular genes amplified 
using primers corresponding to conservative sequences. DNA fragments may 
be cloned preferentially in a piasmid vector (or diluted to the level of one 
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molecule per well in multiwell plates and than amplified in vitro). Clmes 
prepared this way may be resequenced as described above. Two types of 
information are obtained. First of all, a catalogue of different species may be 
defined as wdl as the density of the individuals for each species. Another 
5 segment of informadon may be used to measure the influence of ecological 
factors or pollution on the ecosystem. It may reveal whether some species are 
eradicated or whether the abundance lados among species is altered due to the 
pollution. The method also is applicable for sequencing DNAs from fossils. 

Example 13 

10 DNA Sequencing 

An array of subairays allows for efficient sequencing of a small 
set of samples arrayed in the form of replicated subarrays; For example, 64 
samples may be arrayed on a 8 X-8 mm subanay and 16 X 24 subarrays may 
be replicated on a 15 X 23 cm membrane with 1 mm wide spacers between 

15 the subarrays. Several replica membranes may be made. For example, 
probes from a universal set of three thousand seventy-two 7-mcrs may be 
divided in thirty-two 96-well plates and labelled by Idnasing. Four 
membranes may be processed in parallel during one hybridization cycle. On 
each membrane, 384 probes may be scored. All probes may be scored in two 

20 hybridization cycles. Hybridization intensities may be scored and the 
sequence assembled as described below. 

If a single sample subanay or subarrays contains several 
unknowns, especially when similar samples are used, a smaller number of 
probes may be sufficient if they are intelligently selected on the basis of 

25 results of previously scored probes. For example, if probe AAAAAAA is not 
positive, there is a small chance that any of 8 overlapping probes are positive. 
If AAAAAAA is positive, then two probes are usually positive. The 
sequencing process in this case consists of first hybridizing a subset of 
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mininially overlqyped probes to define positive anchors and then to 
successively select probes which confirms one of the most likely hypotheses 
about the order of anchors and »ze and type of gaps between them. In this 
second i^iase, pools of 2-10 probes may be used where each probe is selected 
S to be positive in only one DNA sample which is diffoent from die samples 
expected to be positive with other probes from the pool. 

The subanay approach aUows efficient implementation of probe 
competition (overlapped probes) or probe cooperation (continuous stacking of 
probes) in solving branching problems. After hybridization of a universal set 
10 of probes the sequence assembly program determines candidate sequence 
subfragments (SFs). For the furtiier assembly of SFs, additional information 
has to be provided (from overlapped sequences of DNA fragments^ similar 
sequences, single pass gel sequences, or from other hybridization or restriction 
mapping data). Competitive hybridization and continuous stacking intmctions 
15 have been proposed for SF assembly. These approaches are of limited 
practical value for sequencing of large numbers of samples by SBH wherein 
a labelled probe is applied to a sample affixed to an array if a uniform array 
is used. Fortunately, analysis of small numbers of samples using replica 
subarrays allows efficient implementation of both approaches. On each of the 
20 replica subarrays, one branching point may be tested for one or more DNA 
samples using pools of probes similarly as in solving mutated sequences in 
different samples spotted in the same subarray (see above). 

If in each of 64 samples described in this example, there are 
about 1(K) branching points, and if 8 samples are analyzed in parallel in each 
25 subarray, then at least 800 subarray probings solve all branches. This means 
Uiat for the 3072 basic probings an additional 800 probings (25%) are 
employed. More preferably, two probings are used for one branching point. 
If the subarrays are smaller, less additional probings are used. For example, 
if subarrays consist of 16 samples, 200 additional probings may be scored 
30 (6%). By using 7-mer probes (N,.2B7N|.2) and competitive or collaborative 
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brancbing solving ai>proaches or both, fragments of about 1000 bp fragments 
may be assembled by about 4000 probings. Furthennoie, using 8-mer prdbcs 
(NB«N) 4 kb or longer fragments may be assembled with 12,000 probings. 
Gapped probes, for example, NB4NB3N or NB4NB4N may be used to reduce 
5 the number of branching points. 

Example 14 

DNA Analysis by Transient Attachment to Subarrays of Pn4>es 
and Ligation of Labelled Probes. 

Oligonucleotide probes having an informative length of four.to 

10 40 bases are synthesized by standard chemistry and stored in tubes or in 
multiweil plates. Specific sets of probes comprising one to 10,000 probes are 
anayed by deposition or in situ synthesis on separate supports or distinci 
sections of a larger support. In the last case, sections or subarrays may be 
separated by physical or hydrophobic barriers. The probe arrays may be 

15 prepared by in situ synthesis. A sample DNA of appropriate size is 
hybridized with one or more specific arrays. Many samples may be 
interrogated as pools at the same subarrays or independently with different 
subarrays within one support. Simultaneously with the sample or subsequently, 
a single labelled probe or a pool of labelled probes is added on each of the 

20 subanays. If attached and labelled probes hybridize back to back on the 
complementary target in the sample DNA they are ligated. Occurrence of 
ligation will be measured by detecting a label from the probe. 

This procedure is a variant of the described DNA analysis 
process in which DNA samples are not permanently attached to the support. 

2S Transient attachment is provided by probes fixed to the support. In this case 
there is no need for a target DNA arraying process. In addition, ligation 
allows detection of longer oligonucleotide sequences by combining short 
labelled probes with short fixed probes. 
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The process has several unique features. Basically, the 
transient attachment of the taiget allows its reuse. After ligation occur the 
target may be released and the label will stay covalently attached to the 
support. This feature allows cycling the target and production of detectable 
S signal with a small quantity of the target. Under typtimal conditions, targets 
do not need to be amplified* e*g. natural sources of the DNA samples may be 
directly used for diagnostics and sequencing purposes. Targets may be 
released by cycling the temperature between efficient hybridization and 
effident melting of duploces. More preferablly, there is no cycling. The 
10 temperature and concentratims of components may be defined to have an 
equilibrium between free targets and targets ent^ed in hybrids at about 
50:50% level. In this case there is a continuous production of ligated 
products. For different purposes different equilibrium ratios are qptiroal. 

An electric field may be used to enhance target use. At the 
15 beginning, a horizontal field pulsing within each subarray may be employed 
to provide for faster target sorting. In this phase, the equilibrium is moved 
toward hybrid formation, and unlabelled probes may be used. After a target 
sorting phase, an appropriate washing (which may be helped by a vertical 
electric field for restricting movement of the samples) may be performed. 
20 Several cycles of discriminative hybrid melting, target harvesting by 
hybridization and ligation and removing of unused targets may be introduced 
to increase specificity. In the next step, labelled probes are added and vertical 
electrical pulses may be applied. By increasing temperature, an optimal free 
and hybridized target ratio may be achieved. The vertical electric field 
25 prevents diffusion of the sorted targets. 

The subarrays of fixed probes and sets of labelled probes 
(specially designed or selected from a universal probe set) may be arranged 
in various ways to allow an efficient and flexible sequencing and diagnostics 
process. For example, if a short fragment (about 100-500 bp) of a bacterial 
30 genome is to be partially or completely sequenced, small arrays of probes (5- 
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30 bases in length) designed on the bases of known sequence may be used. 
If interrogated with a different pool of 10 labdled probes per wbarray, an 
anay of 10 subanays each having 10 probes, allows checking of 200 bases, 
assuming that only two bases connected by ligation are scored. Under the 
5 conditions where mismatches are discriminated throughout the hybrid, probes 
may be displaced by more than one base to cover the longer target with the 
same number of probes. By uang long probes, the target may be interrogated 
directly without amplification or isolation from the rest of DNA in the sample. 
Also, several targets may be analyzed (screened for) in one sample 

10 simultaneously. If the obtained results indicate occurrence of a mutation (or 
a pathogen), additional pools of probes may be used to detect type of die 
mutation or subtype of pathogen. This is a desirable feature of the process 
which may be very cost effective in preventive diagnosis where only a small 
fraction of patients is expected to have an infection or mutation. 

15 In the processes described in the examples* various detection 

methods may be used, for example, radiolabels, fluorescent labels, enzymes 
or antibodies {chemiluminescence), large molecules or particles detectable by 
light scattering or interferometric procedures. 

Example 15 

20 Oligonucleotide Probes and Targets Suitable for SBH 

In order to obtain experimental sequence data defined as a 
matrix of (number of fragments-clones) x (number of probes), the number 
of probes may be reduced depending on the number of fragments used and 
vice versa. The optimal ratio of the two numbers is defined by the 
25 technological requirements of a particular sequencing by hybridization process. 

There are two parameters which influence the choice of probe 
length. The first is the success in obtaining hybridization results that show the 
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requiied degree of discrimination. The second is the tedinological fea^ility 
of synthesis of the required number of probes. 

Hie requirement of obtaining sufficient hybridization 
discrimination with practical and useful amounts of target nucleic acid limits 
5 the probe length. It is difficuh to obtain a sufficient amount of hybrid with 
short probes, and to discriminate end mismatches with long probes. 
Traditionally the use of probes shorter than 11-mers in the Hterature, is 
limited to very stable probes [Estivill et al., Nucl. Acids Res.l5: 1415 (1987)] 
On the other hand, probes longer than IS bases discriminate end mismatches 
10 with difficulty (Wood el al., Proc. Natl. . Acad. Sci. USA 22: 1585 (1985)]. 

One solution for the problems of unstable pr«rt>es and end 
mismatch discrimination is the use of a group of longer probes representing 
a single shorter probe in an informational sense. For example, groups of 
sixteen 10-mers may be used instead of single 8-mers. Every member of the 
15 group has a common core 8-mer and one of three possible variations on outer 
positions with two variations at each end. The probe may be rq)resented as 
5*(A, T, C, G) (A. T, C, G) Bg (A, T, C, G) 3'. With this type of probe one 
does not need to discriminate the non-informative end bases (two on y end, 
and one on 3* end) since only the internal 8-mer is read. This solution 
20 employs a higher mass amounts of probes and label in hybridization reactions. 

These disadvantages are eliminated by the use of a few sets of 
discriminative hybridization conditions for oligomer probes as short as 6-mers. 

The number of hybridization reactions is dependent on the 
number of discrete labelled probes. Therefore in the cases of sequencing 
25 shoner nucleic acids using a smaller number of fragments-clones than the 
number of oligonucleotides, it is better to use oligomers as the target and 
nucleic acid fragment as probes. 

Target nucleic acids which have undefined sequences may be 
produced as a mixture of representative libraries in a phage or plasmid vector 
30 having inserts of genomic fragments of different sizes or in samples prepared 
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by PGR. Inevitable gaps and uncertainties in alignment of sequenced 
fragments arise from nonrandom or rq>etitive sequence organization of 
complex genomes and difficulties in cloning poisonous sequences in 
Escherichia coli. These problems are inherent in sequoidng large complex 
5 molecules using any method. Such problems may be minimized by the choice 
of libraries and number of subclones used for hybridization. Alternatively, 
such difficulties may be ov^come through the use of amplified target 
sequences, e.g. by PGR amplification, ligation reactions, ligation-amplified 
reactions, etc. 

10 Nucleic adds and methods for isolating, cloning and sequencing 

nucleic acids are well known to those of skill in the art. See e.g., Ausube) 
et al,, Current Proiocols in Molecular Biology, Vol. 1-2, John Wiley & Sons 
(1989); and Sambrook et al.. Molecular Cloning A Laboratory Manual, 2nd 
Ed., Vols. 1-3. Cold Springs Harbor Press (1989), both of which are 

15 incorporated by reference herein. 

SBH is a well developed technology that may be practiced by 
a number of methods known to those skilled in the art. Specifically, 
techniques related to sequencing by hybridization of the following documents 
is incorporated by reference herein: Drmanac et al., U.S. Patent No. 

20 5,202,231 (hereby incorporated by reference herein) - Issued April 13, 1993; 
Drmanac et al.. Genomics, 4. 114-128 (1989); Drmanac et al., Proceedings 
of the First Int'L Conf. Electrophoresis Supercomputing Human Genome 
Cantor, DR it Lim HA eds, World Scientific Pub. Co., Singapore, 47-59 
(1991); Drmanac et al.. Science, 260, 1649-1652 (1993); Lehrach et al., 

25 Genome Analysis: Genetic and Physical Mapping, 1, 39-81 (1990), Cold 
Spring Harbor Laboratory Press; Drmanac et al., Nuci Acids Res., 4691 
(1986); Slevanovic et al.. Gene, 79, 139 (1989); Panusku et ai., MoL Biol. 
EvoL, 1, 607 (1990); Nizeiic et al., Nuci Acids Res., 19, 182 (1991); 
Drmanac et al , 7. BiomoL Struct. Dyn. , 5, 1085 (1991); Hoheisel et al.. MoL 

30 Gen., 4, 125-132 (1991); Strezoska etal., Proc: Nat'l. Acad. ScL (USA), 88, 
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10089 (1991); Dnnanac el al., NucL Acids Res., 19. 5839 (1991); and 
Drmanac et al*. Int. J. Genome Res., 1, 59-79 (1992). 

Example 16 

Determining Sequence from Hybridization Data 

Sequence assembly may be inteirupted where ever a given 
overlapping (N-1) mer is duplicated two or more dmes. Then either of the 
two N-mers differing in the last nucleotide may be used in extending the 
sequence. This branching point limits unambiguous assembly of sequence. 

Reassembling the sequence of known oligonucleotides that 
hybridize to the target nucleic add to generate the complete sequence of the 
target nucleic acid may not be accomplished in some cases. This is because 
some information may be lost if the target nucleic acid is not in fragments of 
appropriate size in relation to the size of oligonucleotide that is used for 
hybridizing. The quantity of information lost is proportional to the length of 
a target bdng sequenced. However, if sufficiently short targets are used, their 
sequence msy be unambiguously determined. 

The probable frequency of duplicated sequences that would 
interfere with sequence assembly which is distributed along a certain length 
of DNA may be calculated. This derivation requires the introduction of the 
definition of a parameter having to do with sequence organization: the 
sequence subfragment (SF). A sequence subfragment results if any pan of the 
sequence of a target nucleic acid starts and ends with an (N-l)mer that is 
repeated two or more times within the target sequence. Thus, subfragments 
are sequences generated between two points of branching in the process of 
assembly of the sequences in the method of the invention. The sum of all 
subfragments is longer than the actual target nucleic acid because of 
overlapping short ends. Generally, subfragments may not be assembled in a 
linear order without additional information since they have shared (N-l)mers 
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at their ends and starts. Diffeient numbers of subfragments are obtained for 
each nudeic acid target depending on the number of its repeated (N-1) mers. 
The number depends on the value of N-1 and the length of the taiget. 

Probability calculations can estimate the interrelationship of the 
5 two &ctors. If the ordering of positive N-mers is accomplished by using 
overlapping sequences of loigth N-1 or at an average distance of A^, the N-1 
of a fragment Lf bases long is given by equation one: 
N^=l+Ao xj:kxp(k,l,) 
Where K greater than or = 2» and P (K, L^) represents the probability of an 

10 N-mer occurring K-times on a fragment Lf base long. Also, a computer 
program that is able to form subfragments from the content of N*mers for any 
given sequence is described below in Example 18. 

The number of subfragments increases with the increase of 
lengths of togments for a given length of probe. Obtained subfragments may 

15 not be uniquely ordered among themselves. Although not complete, this 
information is very useful for comparative sequence analysis and the 
recognition of functional sequence characteristics. This type of information 
may be called partial sequence. Another way of obtaining partial sequence is 
the use of only a subset of oligonucleotide probes of a given length. 

20 There may be relatively good agreement between predicted 

sequence according to theory and a computer simulation for a random DNA 
sequence. For instance, for N-l - 7, [using an 8-mer or groups of sixteen 
10-mers of type 5' (A,T,C,G) Bg (A,T,C,G) 3*] a target nucleic acid of 200 
bases will have an average of diree subfragments. However, because of the 

25 dispersion around the mean, a library of target nucleic acid should have 
inserts of SOO bp so that less than 1 in 2000 targets have more than three 
subfragments. Thus, in an ideal case of sequence determination of a long 
nucleic acid of random sequence, a representative library with sufficientiy 
short inserts of target nucleic acid may be used. For such inserts, it is 

30 possible to reconstruct the individual target by the method of the invention. 
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The entire sequofice of a large nucleic add is then obtained by overiapping of 
the defined individual insert sequences. 

To reduce the need for very short fragments, e.g. SO bases for 
8-mer jmbes. The infonnation contained in the overk^ped fragments present 
5 in evwy random DNA liragmentadon process like cloning, or random PGR is 
used. It is also possible to use pools of short physical nucleic acid fragmoits. 
Using 8-mers or ll-mers like 5* (A, T, C, G) N, (A, T. C ,G )3' for 
sequencing 1 megabase, instead of needing 20,000 50 bp fragments only 2, 100 
samples are sufficient. This number consists of 700 random 7 kb clones 
10 (basic library), 1250 pools of 20 clones of 500 bp (subfragmenis ordering 
library) and 150 clones from jumping (or similar) library. The developed 
algorithm (see Example 18) regenerates sequence using hybridization data of 
th these described samples. 

Example 17 

15 Hybridization With Oligonucleotides 

Oligonucleotides were either purchased from Genosys.lnc, 

Houston, Texas or made on an Applied Biosystems 381A DNA synthesizer. 

Most of the probes used were not purified by HPLC or gel electrophoresis. 

For example, probes were designed to have both a single perfectly 
20 complementary target in interferon, a M13 clone containing a 921 bp Eco 

RI-Bgl II human Bl - interferon fragment (Ohno and Tangiuchi, Proc. Natl. 

Acad. Sci. 74: 4370-4374 (1981)]. and at least one target with an end base 

mismatch in M13 vector itself. 

End labelling of oligonucleotides was performed as described 
25 [ManiaUs et al.. Molecular Cloning: A Laboratory Manual, Cold Spring 

Harbor Laboratory Cdld Spring Harbor, New York (1982)] in 10 pd 

containing T4-polynucleotide kinase (5 units Amersham), y*^'^-ATP (3.3 pM, 
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10 ^Ci Amersham 3000 CiymM) and oligonucleotide (4 pM, 10 ng). Specific 
activities of the probes were 2«5-S X 10 9 cpm/nM. 

Single stranded DNA (2 to 4 ^ in 0.5 NaOH» l.S M NaQ) 
was spotted on a Gene Screen membrane wetted with the same solution, the 
S filters were neutralized in O.OS M Na2HP04 pH baked in an oven at 
80^ for 60 min. and UV irradiated for 1 min. Then, the filters woe 
incubated in hybridization solution (0.5 M Na2HP04 pH 7.2» 7% sodium 
lauroyl sarcosine for 5 min at room temperature and placed on the surface of 
a plastic Petri dish. A drop of hybridization solution (10 Ol, 0.5 M Na2HP04 

10 pH 12, 1% sodium lauroy] sarcosine) with a ^P end labelled oligomer probe 
at 4 nM concentration was placed over 1*6 dots per filter, overlaid with a 
square piece of polyethylene (approximately 1 X I cm.), and incubated in a 
moist chamber at the indicated temperatures for 3 hr. Hybridization was 
stopped by placing the filter in 6X SSC washing solution for 3 X 5 minute at 

15 (y*C to remove unhybridized probe. The filter was either dried, or further 
washed for the indicated times and temperatuxes, and autoradiographed. For 
discrimination measurements, the dots were excised from the dried filters after 
autoradiography [a phosphoimager (Molecular Dynamics, Sunnyvale, 
California) may be used) placed in liquid scintillation cocktail and counted. 

20 The uncorrected ratio of c^ms for IF and M13 dots is given as D. 

The conditions reported herein allow hybridizadon with very 
short oligonucleotides but ensure discriminations between matched and 
mismatched oligonucleotides that are complementary to and therefore bind to 
a target nucleic acid. Factors which influence the efficient detection of 

25 hybridization of specific short sequences based on the degree of 
discriminations (D) between a perfectly complementary target and an 
imperfectly complementary target with a single mismatch in the hybrid are 
defined. In experimental tests, dot blot hybridizati(»i of twenty-eight probes 
thai were 6 to 8 nucleotides in length to two M13 clones or to model 
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oligonucleotides bound to membrane filters was accomplished. The principles 
guiding the experimental procedures are given below. 

Oligonucleodde hybridizaticm to filter bound target nucleic 
adds only a few nucleotides longer than the probe in conditioiis of probe 
5 excess is a pseudo-first order reaction with respect to target concentration. 
This reaction is defined by: 

Wherdn S, and are target sequence concentrations at time t and ^ 
respectively. (OP) is probe concentration and t is temperature. The rate 
10 constant for hybrid formation, kfc increases only slightly in the 0^ to 3(fC 
range (Porschke and Eigen, J. Mol, BioL ^2: 361 (1971); Cmg et al., J. 
MoL BioL 62: 383 (1971)], Hybrid melting is a first order reaction with 
respect to hybrid concentration (here replaced by mass due to filter bound 
state) as shown in: 

15 H/H, = e^^* 

In this equation, H, and are hybrid concentrations at times t and to» 
respectively; is a rate constant for hybrid melting which is dependent on 
temperature and salt concentration [Ikuia et al.. Nucl. Acids Res. 15: 797 
(1987); Porsclike and Egen, 7. Mol, BioL §2: 361 (1971); Craig et al., 7. 

.20 MoL BioL 62: 303 (1971)]. During hybridization, which is a strand 
association process, the back, melting, or strand dissociation, reaction takes 
place as well. Thus, the amount of hybrid formed in time is result of forward 
and back reactions. The equilibrium may be moved towards hybrid formation 
by increasing probe concentration and/or decreasing temperature. However, 

25 during washing cycles in large volumes of buffer, the melting reaction is 
dominant and the back reaction hybridization is insignificant, since the pn^e 
is absent. This analysis indicates workable Short Oligonucleotide 
Hybridization (SOU) conditions call be varied for probe concentration or 
temperature. 

30 D or discrimination is defined in equation four: 
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D = H, (U / Hi (O 
Hp (t, ) and Hi (C) are the amounts hybrids remaining after a washing time» 
t^, for the identical amounts of perfectly and imperfectly complementary 
duplex, respectively. For a given temperature, the discriminadoh D changes 
5 with the 10 length of washing time and reaches the maximal value when H^ 
= B which is equation five. 

The background, B, rq>resents the lowest hybridization signal 
detectable in the system. Since any further decrease of H{ may not be 
examined, D increases upon continued wadung. Washing past t» just 
10 decreases relative to B, and is seen as a decrease in D. The optimal 
washing time, t^, for imperfect hybrids, from equation three and equation five 
is: 

= -in (B / Hi (to))/ k., 
Since Hp is being washed for the same , combining equations, one obtains 
15 the optimal discrimination function: - 

J) ^ g h (Bmi(«))km.p/km,i X Hp(to) / B 

The change of D as a function, of T is important because of the choice of an 
optimal washing temperature. It is obtained by substituting the Arhenius 
equation which is: 

into the previous equation to form the final equation: 

D = H,((g/B X (B/H, (g) ' ^ ^-i - ^ V'''' ; 

Wherein B is less tiian H; (Iq). 

Since the activation energy for perfect hybrids, E^p , and the 

25 activation energy for imperfect hybrids, E^j, can be either equal, or E«j less 
than E«p D is temperature independent, or decreases with increasing 
temperature, respectively. This result implies that the search for stringent 
temperature conditions for good discrimination in SOH is unjustified. By 
washing at lower temperatures, one obtains equal or beuer discrimination, but 

30 the time of washing exponentially increases with the decrease of temperature. 
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Discriminadon more strongly decreases with T, if Hjlh) increases relative to 
Hp (to). 

D at Iowa* temperatures ikpends to a higher degree on the Hp 
(tQ)/B ratio than on the Hp (!«) / Hj (Iq) ratio. This result indicates that it is 
5 better to obtain a sufficient quantity of Hp in the hybridization regardless of 
the discrimination that can be achieved in this stq>. Better discrimination can 
then be obtained by washing, since the higher amounts of perfect hybrid allow 
more time for differential melting to show an effect. Similariy, using larger 
amounts of target nucleic add a necessary discrimination can be obtained even 
10 with small differences between K„,p and K^^i. 

Extrapolated to a more complex situation than covered in this 
simple model, the result is that washing at lower temperatures is even more 
important for obtaining discrimination in the case of hybridization of a probe 
having many end*mismatches within a given nucleic acid target. 
IS ' Using the described theoretical principles as a guide for 

experiments, reliable hybridizations have been obtained with prol)es six to 
eight nucleotides in length. AH experiments were performed with a floating 
plastic sheet providing a film of hybridization solution above the filter. -This 
procedure allows maximal reduction in the amount of probe» and thus reduced 
20 label costs in dot blot hybridizations. The high concentration of sodium 
lauroyl sarcosine instead of sodium lauroyl sulfate in the phosphate 
hybridization buffer allows dropping the reaction from room temperature 
down to ^''C. Similarly, the 4-6 X SSC, 10% sodium lauroyl sarcosine 
buffer allows hybridization at temperatures as low as The detergent in 
25 these buffers is for obtaining tolerable background with up to 40 nM 
concentrations of labelled probe. Preliminary characterization of the thermal 
stability of short oligonucleotide hybrids was determined on a prototype 
octamer with 50% G+C contoit, i.e. probe of sequence TGCTCATG. The 
theoretical expectation is that this probe is among the less stable octamers. 
30 Its transition enthalpy is similar to those of more stable heptamers or, even to 
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probes 6 nucleotides in length (Bressiauer et al., Proc. Natl Acad. ScL 
V.S.A. aj: 3746 (1986)). Parameter T^, the temperature at which 50% of the 
hybrid is melted in unit time of a minute is IS'C. Hie result shows that 
is \y<l lower for the 8 bp hybrid than for an 11 bp duplex [Wallace et al., 
5 Nud^c Adds Res. fi: 3543 (1979)]. 

In addition to experiments with model oligonucleotides, an M13 
vector was chosen as a system for a practical demonstration of short 
oligonucleotide hybridization. The main aim was to show useful 
end-mismatch discrimination with a target similar to the ones which will be 

10 used in various apjdications of the method of the invendon. Oligonucleotide 
probes for the M13 model were chosen in such a way that the M13 vector 
itself contains the end mismatched base. Vector IF, an M13 recombinant 
containing a 921 bp human interferon gene insert, carries single perfectly 
matched target. Thus, IF has either the identical or a higher number of 

15 mismatched targets in comparison to the M13 vector itself. 

Using low temperature conditions and dot blots, sufficient 
differences in hybridization signals were obtained between lie dot containing 
the perfect and the mismatched targets and the dot containing the mismatched 
targets only. This was true for the 6-mer oligonucleotides and was also true 

20 for the 7 and 8-mer oligonucleotides hybridized to the large IF-M13 pair of 
nucleic acids. 

The hybridization signal depends on the amount of target 
available on the filter for reaction with the probe. A necessary control is to 
show that the difference in sign intensity is not a reflection of varying amounts 

25 of nucleic acid in the two dots. Hybridization with a probe that has the same 
number and kind of targets in both IF and Ml 3 shows that there is an equal 
amount of DNA in the dots. Since the efficiency of hybrid formation 
increases with hybrid length, the signal for a duplex having six nucleotides 
was best detected with a high mass of oligonucleotide target bound to the 

30 filter. Due to their lower molecular weight, a larger number of 
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oligonucleotide taiget molecules can be bound to a given surface area wben 
compared to huge molecules of nucleic acid that serves as target. 

To measure the sensitivity of detection with unpurified DNA, 
various amounts of phage siq)ematants were spotted on the filter and 
5 hybridized with a ^P-labelled ociamer. As litde as SO million unpurified 
phage containing no more than 0.5 ng of DNA gave a detectable rignal 
indicating that sensitivity of the short oligonucleotide hybridization method is 
sufficient. Reacdon time is short, adding to the practicality. 

As mentioned in the theoretical section above, the equilibrium 
10 yield of hybrid depends oil probe concentration and/or temperature of 
reaction. For instance, the signal level for the same amount of taiget with 4 
nM octamer at IS'^C is 3 times lower than with a probe concentration of 40 
nM, and is decreased 4.S-times by raising the hybridization temperature to 
25^. 

15 The utility of the low temperature wash for achieving maximal 

discrimination is demonstrated. To make the phenomenon visually obvious, 
50 times more DNA was put in the M13 dot than in the IF dot using 
hybridization with a vector specific probe. In this way, the signal after the 
hybridization step with the actual probe was made stronger in the, mismatched 
20 that in the matched case. The Hp /Hj ratio was 1:4. Inversion of signal 
intensities after prolonged washing at TC was achieved without a massive loss 
of perfect hybrid, resulting in a ratio of 2: 1. In contrast, ii is impossible to 
achieve any discrimination at 25''C, since the matched target signal is already 
brought down to the background level with 2 minute washing; at the same 
25 time, the signal from the mismatched hybrid is still detectable. The loss of 
discrimination at 13°C compared to 7*C is not so great but is clearly visible. 
If one considers the 90 minute point at TQ and the 15 minute point at 13*C 
when, the mismatched hybrid signal is near the background level, which 
repres^ts optimal washing times for the req^ective conditions, it is obvious 
30 that the amount of several times greater at TC than at 13**C. To illustrate this 
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fiirther, the time course of the change discrimination with washing of the same 
amount of starting hybrid at the two temperatures shows the higher maximal 
D at tiie lower temperature. These results confirm the trend in the change of 
D with tenq)eratuie and the ratio of amounts of the two types of hybrid at the 
5 start of the washing step. 

In order to ^ow the geneial utility of the short oligonucleotide 
hybridization conditions, we have looked hybridization of 4 hqitamers, 10 
octamers and an additional 14 probes up to 12 nucleotides in length in our 
simple M13 system. These include-the nonamer GnTTTTAA and octamer 

10 GGCAGGCG representing the two extremes of GC content. Although GC 
content and sequence are expected to influence the stability of short hybrids 
[Kesslauer et al., Proc. Natl Acad. ScL U.S.A. fi2: 3746 (1986)], the low 
temperature short oligonucleotide conditions were applicable to all tested 
probes in achieving sufficient discrimination. Since the best discrimination 

15 value obtained with probes 13 nucleotides in length was 20, a several fold 
drop due to sequence variation is easily tolerated. 

The Ml 3 system has the advantage of showing the effects of 
target DNA complexity on the levels of discrimination. For two octamers 
having either none or five mismatched targets and differing in only one GC 

20 pair the observed discriminations were 18.3 and 1.7, respectively. 

In order lo show the utility of this method, three probes 8 
nucleotides in length were tested on a collection of 51 plasmid DNA dots 
made from a library in Bluescript vector. One probe was present and specific 
for Bluescript vector but was absent in Ml 3, while the other two probes had 

25 targets that were inserts of known sequence. This system allowed the use of 
hybridization negadve or positive control DNAs with each probe. This probe 
sequence (CTCCCTTT) also had a complementary target in the interferon 
insert. Since the M13 dot is negative while the interferon insert in either Ml 3 
or Bluescript was positive, the hybridization is sequence specific. Similarly, 

30 probes that detect the target sequence in only one of 51 inserts, or in none of 
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Ihe examined inserts along with controls that confirm that hybridization would 
have occurred if the appropriate targets were present in the clones. 

Thermal stability curves for very short oligonucleotide hybrids 
that are 6-8 nucleotides in length are at least 15*YI low^ than for hybrids 
S 1 1-12 nucleotides in length [Fig. 1 and Wallace et al., Nucleic Adds Res. 6: 
3543*3557 (1979)]. However^ performing the hybridization reaction at a low 
temperature and with a vory practical 0.4-40 nM concentration of 
oligonucleotide probe allows the detection of complementary sequence in a 
known or unknown nucleic add target. To determine an unknown nucleic 
10 add sequence completdy, an entire set containing 65»53S 8-mer probes may 
be used. Sufficient amounts of nucleic acid for this purpose are present in 
convenient biological samples such as a few microliters of Ml 3 culture, a 
plasmid prep from 10 ml of bacterial culture or a single colony of bacteria, 
or less than 1 ^1 of a standard PCR reaction. 
15 Short oligonucleotides 6-10 nucleotides long give excellent 

discrimination. The relative decrease in hybrid stability with a single end 
mismatch is greater than for longer probes. Results with the octamer 
TGCTCATG support this conclusion. In the experiments, the target with a 
G/T end mismatch, hybridization to the target of this type of mismatch is the 
20 most stable of all other types of oligonucleotide. This discrimination achieved 
is the same as or greater than an internal G/T mismatch in a 19 base paired 
duplex greater than an internal G/T mismatch in a 19 paired duplex flkuta et 
al., NucL Acids Res. 15: 797 (1987)]. Exploiting these discrimination 
properties using the described hybridization conditions for short 
25 oligonucleotide hybridization allows a very precise determination of 
oligonucleotide targets. 

In contrast to the ease of detecting discrimination between 
perfect and imperfect hybrids, a problem that may exist with using very short 
oligonucleotides is the prqiaration of sufficient amounts of hybrids. In 
30 practice, the need to discriminate Hp and is aided by increasing the amount 
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of DNA in the dot and/or the probe conoentzationy or by decreasing the 
hybiidixation tmq)erature. However, higher probe concentrations usually 
increase background. Moreover, there are limits to the amounts of target 
nucldc add that are practical to use. This problems was solved by the higher 
S concentration of the detergent Sarcosyl whidi gave an effective background 
with 4 nM of probe. Further improvements may be effected dther in the use 
of competitors for unspedfic binding of probe to filter, or by changing the 
hybridization siq>port material. Moreover, for probes having less than 45 
Kcal/mol (e.g. for many heptamers and a minority of hexamers, modified 

10 oligonucleotides give a more stable hybrid TAsseline, et al. , Proc. Nat 7 Acad. 
ScL £1: 3297 (1984)] dian tiidr unmodified counterparts. The hybridization 
conditions described in this invention for short oligonucleotide hybridization 
using low temperatures give better discriminating for all sequences and duplex 
hybrid inputs. The only price paid in achieving uniformicy in hybridization 

15 conditions for different sequences is an increase in washing time from minutes 
to up to 24 hours depending on the sequence. Moreover, the washing time 
can be further reduced by decreasing the salt concentration. 

Although there is excellent discrimination of one matched 
hybrid over a mismatched hybrids, in short oligonucleotide hybridization, 

20 signals from mismatched hybrids exist, with the majority of the mismatch 
hybrids resulting from end mismatch. This may limit insert sizes that may be 
effectively examined by a probe of a certain length. 

The influence of sequence complexity on discrimination cannot 
be ignored. However, the complexity effects are more significant when 

25 defining sequoice information by short oligonucleotide hybridization for 
specific, nonrandom sequences, and can be overcome by using an appropriate 
probe to target length ratio. The lengtii ratio is chosen to make unlikely, on 
statistical grounds, the occurrence of specific sequences which have a number 
of end-mismatdies which would be able to diminate or falsely invert 

30 discrimination. Results suggest the use of oligonucleotides 6, 7, and 8 
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nucleotides in length on target nucleic add inserts shorter than 0.6, and 
10 kb, respectively. 

Example 18 

Sequendng a Target Using Octamers and Nonamers 

In this example, hybridization conditions that were used are 
described supra in Example 17. Data resulting from the hybridizaticm of 
octamer and nonamer oligonucleotides shows that sequencing by hybridization 
provides an extremely high degree of accuracy. In this experiment, a known 
sequence was used to predict a series of contiguous overiapping component 
octamer and nonamer oligonucleotides. 

In addition to the perfectly matching oligonucleotides, mismatch 
oligonucleotides, mismatch oligonucleotides wherein internal or end 
mismatches occur in the duplex formed by die oligonucleotide and the target 
were examined. In these analyses, the lowest practical temperature was used 
to maximize hybridization formadon. Washes were accomplished at the same 
or lower temperatures to ensure maximal discrimination by utilizing the 
greater dissociation rate of mismatch versus matched oligonucleotide/target 
hybridization. These conditions are shown to be applicable to all sequences 
although the absolute hybridization yield is shown to be sequence dependent. 

The least destabilizing mismatch that can be postulated is a 
simple end mismatch, so that the test of sequencing by hybridization is the 
ability to discriminate perfectly matched oligonucleotide/target duplexes from 
end-mismatched oligonucleodde/target duplexes. 

The discriminative values for 102 of lOS hybridizing 
oligonucleotides in a dot blot format were greater than 2 allowing a highly 
accurate generation of the sequence. This system also allowed an analysis of 
the effect of sequence on hybridization formation and hybridization instability. 
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Qne hundred base pairs of a known poition of a human 
0-mVsrfcnm genes prepared by PGR, i.e. a 100 bp target sequence, was 
generated with data resulting from the hybridization of lOS oligonucleotides 
probes of known sequence to the target nucldc acid. The oligonucleotide 
5 probes used included 72 octamer and 21 nonamer oligonucleotides whose 
sequence was perfectly complementary to the target. The set of 93 probes 
provided consecutive overlapping frames of the target sequence e displaced by 
one or two bases. 

To evaluate the effect of mismatches, hybridization was 

10 examined for 12 additional probes that contained at least one end mismatch 
when hybridized to the 100 bp test target sequence. Also tested was the 
hybridization of twelve jmbes with target end-mismatched to four other 
control nucleic acid sequences chosen so that the 12 oligonucleotides formed 
perfectly matched duplex hybrids with the four control DNAs. Thus, the 

15. hybridization of internal mismatched, end-mismatched and perfectly matched 
duplex pairs of oligonucleotide and target were evaluated for each 
oligonucleotide used in the experiment. The effect of absolute DNA target 
concentration on the hybridization with the test octamer and nonamer 
oligonucleotides was determined by defining target DNA concentration by 

20 detecting hybridization of a different oligonucleotide probe to a single 
occurrence non- target site within the co-amplified plasmid DNA. 

The results of this experiment showed that all oUgonucleotides 
containing perfect matching complementary sequence to the target or control 
DNA hybridized more strongly than those oligonucleotides having 

25 mismatches. To come to this conclusion, we examined Hp and D values for 
each probe. defines the amount of hybrid duplex formed between a test 
target and an oligonucleotide probe. By assigning values of between 0 and 10 
to the hybridization obtained for the 105 probes, it was apparent that 68.5% 
of the 105 probes had an greater than 2. 
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Discrimination (D) values were obtained where D was defined 
as the ratio of signal intensities between 1) the dot containing a perfect 
matched diq^lex formed between test oligonucleotide and target or control 
nucleic acid and 2) the dot containing a mismatch duplex formed between the 
5 same oligonucleotide and a different site within the target or control nucleic 
add. Variations in the value of D result from either 1) perturbations in the 
hybridization efficiency which allows visualization of signal over background, 
or 2) the type of mismatch found between the test oligonucleotide and the 
target. The D values obtained in this experimoit were between 2 and 40 for 
10 102 of the 105 oligonucleotide probes examined. Calculations of D for the 
group of 102 oligonucleotides as a whole showed the average D was 10.6. 

There were 20 cases where oligonucleotide/target duplies 
exhibited an end^mismatch. In five of these, D was greater than 10. The 
large D value in these cases is most likely due to hybridization destabiUzatton 
15 caused by other than the most stable (G/T and G/A) end mismatches. The 
other possibility is there was an error in the sequence of either the 
oligonucleotides or the target. 

Error in the target for probes with low Hp was excluded as a " 
possibility because such an error would have affected the hybridization of each 
20 of the other eight overlapping oligonucleotides. There was no apparent - 
instability due to sequence mismatch for the other overlapping 
oligonucleotides, indicating the target sequence was correct. Error in the 
oligonucleotide sequence was excluded as a possibility after the hybridization 
of seven newly synthesized oligonucleotides was re-examined. Only 1 of the 
25 seven oligonucleotides resulted in a better D value. Low hybrid formation 
values may result from hybrid instability or from an inability to form hybrid 
duplex. An inability to form hybrid duplexes would result from either 1) self 
complementarity of the chosen probe or 2) target/target self hybridization. 
Oligonucleotide/oligonucleotide duplex formation may be favored over 
30 oligonucleotide/target hybrid duplex formation if the probe was 
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self-complementary. SimUarly , target/isrgei association may be favored if the 
target was sdf-coraplemeniary or may form internal palindromes. In 
evaluating these possibilities* it was a^Jparent from probe analyus that the 
questionable probes did not form hybrids with themselves. Moreov^, in 
5 examining the contribution of target/target hybridization, it was determined 
that one of tiie questionable oligonucleotide probes hybridized ineffidentiy 
witii two different DNAs containing tiie same target. The low probabiUty that 
two different DNAs have a self-complementary region for the same target 
sequence leads to die conclusion that target/target hybridization did not 

10 contribute to low hybridization formation. Thus, tiiese results indicate tiiat 
hybrid instability and not die inability to form hybrids was die cause of the 
low hybrid formation observed for specific oligonucleotides. The results also 
indicate tiiat low hybrid formation is due to die specific sequences of certain 
oligonucleotides. Moreover, die results indicate diat reliable results may be 

IS obtained to generate siequoices if octamer and nonamer oligonucleotides are 
used. 

These results show thai using the methods described long 
sequences of any specific target nucleic acid may be generated by maximal 
and unique overlap of consuiueni oligonucleotides. Such sequencing methods 

20 are dependent on the content of the individual component oligomers regardless 
of their frequency and their position. 

The sequence which is generated using the algorithm described 
below is of high fidelity. The algorithm tolerates false positive signals from 
the hybridization dots as is indicated from Uie fact die sequence generated 

25 from die 105 hybridization values, which included four less reliable values, 
was correct. This fidelity in sequencing by hybridization is due to the "all or 
none" kinetics of short oligonucleotide hybridization and die difference in 
duplex stability that exists between perfectly matched duplexes and 
mismatched duplexes. The ratio of duplex stability of matched and 

30 end-mismatched duplexes increases witii decreasing duplex length. Moreover. 
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binding energy decieases with decreasing duplex length resulting in a lower 
hybridization efficiency. However, the results provided show that octamer 
hybridization allows the balancing of the factors affecting duplex stability and 
discrinunation to produce a highly accurate method of sequencing by 
5 hybridization. Results presented in other examples show that oligonucleotides 
that are 6, 7, or 8 nucleotides can be effectively used to generate reliable 
sequence on targets that are 0.5 kb (for hexamers) 2 kb (for septamers) and 
6kb (for octamers). The sequoice of long fragments may be overlapped to 
generate a complete genome sequence. 
10 An algorithm to determine sequence by hybridization is 

described in Example 18. 

Example 19 
Algorithm 

This example describes an algorithm for generation of a long 
15 sequence written in a four letter alphabet from constituent k-tupie words in a 
minimal number of sqparate, randomly defined fragments of a starting nucleic 
acid sequence where K is the length of an oligonucleotide probe. The 
algorithm is primarily intended for use in the sequencing by hybridization 
(SBH) process. The algorithm is based on subfragments (SF), infonnative 
20 fragments (IF) and the possibility of using pools of physical nucleic sequences 
for defining informative fragments. 

As described, subfragments may be caused by branch points in 
the assembly process resulting from the repetition of a K-1 oligomer sequence 
in a target nucleic acid. Subfragments are sequence fragments found between 
25 any two rq[}etitive words of the length K-1 that occur in a sequence. Multiple 
occurrences of K-1 words are the cause of interruption of ordering the overlap 
of K- words in the process of sequence generation. Interruption leads to a 
sequence remaining in the form of subfragments. Thus, the unambiguous 
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segmoits between blanching points whose order is not uniquely determined 
are called sequoice subfragments. 

Informative fragments are defined as fragments of a sequence 
that are determined by the nearest ends of overlapped physical sequence 
S fragments. 

A certain number of physical fragments may be pooled without 
losing the possibility of defining informative fragments. The total length of 
randomly pooled fragments depends on the length of k-t\q>les that are used in 
the sequencing process. 

10 The algorithm consists of two main units. The first part is used 

for generation of subfragments from the set of k-tuples contained in a 
sequence. Subfragments may be generated within the coding region of 
physical nucleic acid sequence of certain sizes » or within the informative 
fragments defined within long nucleic acid sequences. Both types of 

IS fragments are members of the basic library. ' This algorithm does not describe 
the determination of the content of the k-tuples of the informative fragments 
of the basic library, i.e. the step of preparation of informative fragments to 
be used in the sequence generation process. 

The second part of the algorithm determines the linear order of 

20 obtained subfragments with the purpose of regenerating the complete sequence 
of the nucleic acid fragments of the basic library. For this purpose a second, 
ordering library is used, made of randomly pooled fragments of the starting 
sequence. The algorithm does not include the step of combining sequences 
of basic fragments to regenerate an entire, megabase plus sequence. This may 

25 be accomplished using the link-up of fragments of the basic library which is 
a prerequisite for informative fragment generation. Alternatively, it may be 
accomplished after generation of sequences of fragments of the basic library 
by this algorithm, using search for their overlap, based on the presence of 
common end-sequences. 
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The algorithm requires neither knowledge of the numb^ of 
sqypeanuices of a given k-tuple in a nucleic acid sequence of the basic and 
ofdoing libraries, nor does it require the information of which k-tuple words 
are presnt on the ends of a fragmmt. The algorithm qpeiaies with the mixed 
5 content of k-tuples of various length. The concept of the algorithm enables 
operations with the k-tuple sets that contain &lse positive and false negative 
k- tuples. Only in specific cases does the content of the false k-tuples 
primarily influence the completeness and correctness of the generated 
sequence. The algorithm may be used for optimization of param^ers in 
10 simulation experiments, as well as for sequence generation in the actual SBH 
experiments e.g. generation of the genomic DNA sequence. In optimization 
of parameters, the choice of the oligonucleotide probes (k*tuples) for practical 
and convenient fragments and/or the choice of the optimal lengths and the 
number of fragments for the defined probes are especially important. 

This part of the algorithm has a central role in the process of 
the generation of the sequence from the content of k-tuples. It is based on the 
unique ordering of k-tuples by means of maximal overlap. The main obstacles 
in sequence generation are specific repeated sequences and false positive * 
and/or negative k-iuples. The aim of this pan of the algorithm is to obtain tiie 
20 minimal number of the longest possible subfragments, with correct sequence. 
This part of the algorithm consists of one basic, and several control steps. A 
two-stage process is necessary since certain information can be used only after 
generation of all primary subfragments. 

The main problem of sequence generation is obtaining a 
25 repeated sequence from word contents that by definition do not carry 
information on the number of occurrences of the particular k-tuples. The 
concept of the entire algorithm depends on the basis on which this problem is 
solved. In principle, Uiere are two opposite approaches: 1) repeated 
sequences may be obtained at the beginning, in the process of generation of 
30 pSFs, or 2) repealed sequences can be obtained later, in the process of the 
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final ordering of the subfragments. In the first case» pSFs contain an excess 
of sequmces and in the second case, they contain a deficit of sequences. The 
first approach requires dimination of the excess sequences generated, and the 
second requires permitting multiple use of some of the subfragments in the 
5 process of the final assembling of the sequence. 

The difference in the two aj^roaches in the degree of strictness 
of the rule of unique overlap of k*tuples. The less severe rule is: k-tuple X 
is unambiguously maximally overlapped with k-tuple Y if and only if, the 
rightmost k-1 end of k-tuple X is present only on the leftmost end of k-tuple 

10 Y. This rule allows the generation of repetitive sequences and the formation 
of suiplus sequences. 

A stricter rule which is used in the second approach has an 
addition caveat: k*tuple X is unambiguously maximally overl^iped with 
k^tupIe Y if and only if, the rightmost K-l end of k-tuple X is present only on 

15 the leftmost end of k-tuple Y and if the leftmost K-l end of k- tuple Y is not 
present on the rightmost end of any other k-tuple. The algorithm based on the 
stricter rule is simpler, and is described herein. 

The process of elongation of a given subfragment is stopped 
when the right k-l end of the last k-tuple included is not present on the left 

20 end of any k-tuple or is present on two or more k-luples. If it is present on 
only one k-tuple the second pan of the rule is tested. If in addition there is 
a k-tuple which differs from the previously included one, the assembly of the 
given subfragment is terminated only on the first leftmost position. If this 
additional k-tuple does not exist, the conditions are met for unique k-l overlap 

25 and a given subfragment is extended to the right by one element. 

Beside the basic rule, a supplementary one is used to allow the 
usage of k-tuples of different lengths. The maximal overlap is the length of 
k-l of the shorter k-tuple of the overlapping pair. Generation of the pSFs is 
performed starting from the first k-tuple from the file in which k-tuples are 

30 displayed randomly and independently from their order in a nucleic acid 
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sequence. Thus, the first k-tuple in the file is not necessarily on the beginning 
of the sequence, nor on the start of the particular subfragment. The process 
of subftagment generation is performed by onteing the k*tuples by means of 
unique ov^lap, which is defined by the described rule. Each used k-tuple is 
erased fiom the file. At the point when there are no further k-tuples 
unambiguously overlapi»ng with the last one included, the building of 
subiragment is tmninated and the buildup of another pSF is started. Since 
generation of a majority of subfragments does not begin from thdr actual 
starts, the formed pSF are added to the k-tuple file and are considmd as a 
longer k-tuple. Another possibility is to form subfragments going in both 
directions from the starting k- tuple. The process ends when further overlap, 
i.e. the extension of any of the subfragments, is not possible. 

The pSFs can be divided in three groups: 1) Subfragments of 
the maximal length and correct sequence in cases of exact k-tuple set; 2) 
short subfragments, formed due to the used of the maximal and unambiguous - 
overlap rule on the incomplete set, and/or the set with some false positive 
k-tuples; and 3) pSFs of an incorrect sequence. The incompleteness of the 
set in 2) is caused by false negative results of a hybridization experiment, as 
weJl as by using an incorrect set of k-tuples. These are formed due to the 
false positive and false negative k-tuples and can be : a) misconnecied 
subfragments; b) subfragments with the wrong end; and c) false positive 
k-tuples which appears as false minimal subfragments. 

Considering false positive k-tuples. there is the possibility for 
the presence of a k- tuple containing more than one wrong base or containing 
one wrong base somewhere in the middle, as well as the possibility for a 
k-tuple with a wrong base on the end. Generation of short, erroneous or 
misconnected subfragments is caused by the latter k-tuples. The k-tuples of 
the former two kinds represent wrong pSFs with length equal to k-tuple 
length. 
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In the case of one folse negative k-tuple, pSFs axe generated 
because of the impossibility of maximal overlapping. In the case of the 
presence of one false positive k-tuple with the wrong base on its leftmost or 
rightmost end, pSFs are generated because of the impossibility of 
5 unambiguous overlapping. When both £alse positive and false negative 
k-tuples with a common k-1 sequence are present in the file, pSFs are 
generated, and one of these pSFs contains the wrong k-tuple at the relevant 
end. 

The process of correcting subfragments with errors in sequence 

10 and the linking of unambiguously connected pSF is performed after 
subfragment generation and in the process of subfragment ordering. The first 
step which consists of cutting the misconnected pSFs and obtaining the final 
subfragments by unambiguous connection of pSFs is described below. 

There are two approaches for the formation of misconnected 

IS subfragments. In the first a mistake occurs when an erroneous k-tuple appears 
on the points of assembly of the repeated sequences of lengths k-1. In the 
second, the repeated sequences are shorter than k-1. These situations can 
occur in two variants each. In the first variant, one of the repeated sequences 
represents the end of a fragment. In the second variant, the repeated sequence 

20 occurs at any position within the fragment. For the first possibility, the 
absence of some k-tuples from the file (false negatives) is required to 
generate a misconnection. The second possibility requires the presence of 
both false negative and false positive k-tuples in the file. Considering the 
repetitions of k-1 sequence, the lack of only one k-tuple is sufficient when 

25 either end is repeated internally. The lack of two is needed for strictiy 
internal repetition. The reason is that the end of a sequence can be considered 
informatically as an endless linear array of false negative k-tupies. From the 
"smaller than k-1 case", only the repeated sequence of the length of k-2, 
which requires two or three specific erroneous k-tuples, will be considered. 
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It is very likely that these will be the only cases which will be detected in a 
real ^perinient, the others being much less frequent. 

Recognition of the misconnected subfragments is more strictly 
defined whm a repeated sequence does not appear at the end of the fragment. 
S In this situation, one can detect further two subiiagments, one of which 
contains on its leftmost, and the other on its rightmost end k-2 sequences 
which are also present in the misconnected subfragment. When the rqieated 
sequence is on the end of the fragment, there is only one subfragment which 
contains k-2 sequence causing the mistake in subfragment formation on its 
10 leftmost or rightmost end. 

The removal of misconnected subframents by their cutting is 
performed according to the common rule: If the leftmost or rightmost 
sequence of the length of k-2 of any subfragments is present in any other 
subfragment, the subfragment is to be cut into two subfragments, each of them 
containing k-2 sequence. This rule does not cover rarer situations of a 
repeated end when there are more than one false negative k-tuple on the point 
of repeated k-1 sequence. Misconnected subfragments of this kind can be 
recognized by using the information from the overlapped fragments, or 
informative fragments of both the basic and ordering libraries. In addition, 
the misconnected subfragment will remain when two or more false negative 
k*tuples occur on both positions which contain the identical k-I sequence. 
This is a very rare situation since it requires at least 4 specific false k-tuples. 
An additional rule can be introduced to cut these subfragments on sequences 
of length k if the given sequence can be obtained by combination of sequences 
shorter than k-2 from the end of one subfragment and the start of another. 

By strict application of the described rule, some completeness 
is lost to ensure the accuracy of the output. Some of the subfragments will 
be cut although they are not misconnected since they fit into the pattern of a 
misconnected subfragment. There are several situations of this kind. For 
example, a fragment, beside at least two identical k*l sequences, contains any 
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k*2 sequence from k-1 or a fragment contains k-2 sequence repeated at least 
twice and at least one false negative k-tuple containing given k-2 sequence in 
the middle, etc. 

The aim of this part of the algorithm is to reduce the number 
S of pSFs to a minimal number of longer subfragments with correct sequence. 
The generation of unique longer subfragments or a complete sequence is 
possible in two situations. The first situation concerns the specific order of 
repeated k-I words. There are cases in which some or all maximally extended 
pSFs (the first group of pSFs) can be uniquely ordered. For example, in 

10 fragment S-Rl-a-R2-b-Rl-c-R2-E where S and E are the start and end of a 
fragment, a, b , and c are different sequences speciftc to respective 
subfragments and Rl and R2 are two k-1 sequences that are tandemly 
repeated, five subfragments are generated (S-Rl . Rl-a-R2. R2-b-Rl , R1-C-R2, 
and R-E). They may be ordered in two ways; the original sequence above or 

15 S-Rl- c-R-b-Rl-a-R-E. In contrast, in a fragment with the same number and 
types of repeated sequences but ordered differently, i.e. S-Rl-a-Rl-b-R-c-R-E, 
there is no other sequence which includes all subfragments. Examples of this 
type can be recognized only after the process of generation of pSFs. They 
represent the necessity for two steps in the process of pSF generation. The 

20 second situation of generation of false short subfragments on positions of 
nonrepeated k-1 sequences when the files contain false negauve and/or 
positive k-tuples is more important. 

The solution for both pSF groups consists of two parts. First, 
the false positive k- tuples appearing as the nonexisting minimal subfragments 

25 are eliminated. All k-iuple subfragments of length k which do not have an 
overlap on either end, of the length of longer than k-a on one end and longer 
than k-b on the other end, are eliminated to enable formation of the maximal 
number of connections. In our experiments, the values for a and b of 2 and 
3, respectively, appeared to be adequate to eliminate a sufficient number of 

30 false positive k-tuples. 
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The meiging of subfiagments that can be uniquely connected 
is accomplished in the second step. The rule for connection is: two 
subfiagments may be unambiguously connected if, and only if, the 
overlapping sequence at the relevant end or start of two subfragments is not 
present at the start and/or end of any other subfragment. 

The exoqjtion is if one subfragment from the considered pair 
has die identical beginning and end. In that case connection is permitted, even 
if there is another subfragment with the same end presmt in the file. The 
main problem here is the precise defmitiDn of overlapping sequence. The 
connection is not permitted if the overlapping sequence unique for only one 
pair of subfragments is shorter than k-2, of it is k-2 or longer but an 
additional subfragment exists with the overlapping sequence of any length 
longer than k-4. Also, both the canonical ends of pSFs and the ends after 
omitting one (or few) last bases are considered as the overlapping sequences. 

After this step some false positive k-tuples (as minimal 
subfragments) and some subfragments with a wrong end may survive. In 
addition, in very rare occasions where a certain number of some specific false 
k-iuples are simultaneously present, an erroneous connection may take place. 
These cases will be detected and solved in the subfragment ordering process, 
and in the additional control steps along with the handling of uncut 
"misconnected" subfragments. 

The short subfragments that are obtained are of two kinds. In 
the common case» these subfragments may be unambiguously connected 
among themselves because of the distribution of repeated k-l sequences. This 
may be done after the process of generation of pSFs and is a good example 
of the necessity for two steps in the process of pSF generation. In the case 
of using the file containing false positive and/or false negative k-tuples, short 
pSFs are obtained on the sites of nonrepeated k-l sequences. Considering 
false positive k-tuples, a k-tuple may contain more than one wrong base (or 
containing one wrong base somewhere in the middle), as well as k-tuple on 
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the end. Generation of short and erroneous (or misconnected) subfragments 
is caused by the latter k-tuples. The k-tuples of the former kind represent 
wrong pSFs with imgth equal to k- Vtxplt length. 

The aim of m^ ging pSF part of the algorithm is the reduction 
5 of the number of pSFs to the minimal number of longer subfragments with the 
correct sequence. All k-tuple subfragments that do not have an overlap on 
either end, of the length of longer than k-a on one» and longer than k-b on the 
other end, are eliminated to enable the maximal number of connections. In 
this way, the majority of false positive k-tuples are discarded. The rule for 

10 connection is: two subfragments can be unambiguously connected if, and only 
if the overlapping sequence of the relevant end or start of two subfragments 
is not present on the start and/or end of any other subftagment. The 
exception is a subfragment with the identical beginning and end. In that case 
connection is permitted, provided that there is another subfragment with the 

IS same end present in the file. The main problem here is of precise definition 
of overlapping sequence. The presence of at least two specific false negative 
k- tuples on the points of repetition of k-1 or k-2 sequences, as well as 
combining of the false positive and false negative k-tuples may destroy or 
"mask" some overlapping sequences and can produce an unambiguous, but 

20 wrong connection of pSFs. To prevent this, completeness must be sacrificed 
on account of exactness: the connection is not permitted on the end-sequences 
shorter than k-2, and in Uie presence of an extra overlapping sequence longer 
than k-4. The overlapping sequences are defined from the end of the pSFs, 
or omitting one, or few last bases. 

25 In the very rare situations, with the presence of a certain 

number of some specific false positive and false negative k-tuples, some 
subfragments with the wrong end can survive, some false positive k-tuples (as 
minimal subfragments) can remain, or the erroneous connection can take 
place. These cases are detected and solved in the subfragments ordering 
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process, and in the additional control steps along with the handling of uncut, 
misconnected subiragments. 

The process of oid^g of subfragments is similar to die 
process of their generation. If one considers subfragments as longer k-tuples, 
5 ordering is performed by their unambiguous connection via overlapping ends. 
The informational basis for unambiguous connection is the division of 
subfragments generated in fragments of the basic library into groups 
npiesrating segments of those fragments. The method is analogous to the 
biochemical solution of this problem based on hybridization with longer 
10 oligonucleotides with relevant connecting sequence. The connecting sequences 
are generated as subfragments using the k-tupie sets of the appropriate 
segments of basic library fragments. Relevant segments are defined by the 
fragments of the ordering library that overlap with the respective fragments 
of the basic library^ The shortest segments are informative fragments of the 
15 ordering library. The longer ones are several neighboring informative 
fragments or total overlapping portions of fragments corresponding of the 
ordering and basic libraries. In order to decrease the number of separate 
samples, fragments of the ordering library are randomly pooled, and the 
unique k- tuple content is determined. 
20 By using the large number of fragments in the ordering library 

very shon segments are generated, thus reducing the chance of the multiple 
appearance of the k-1 sequences which are the reasons for generation of the 
subfragments. Furthermore, longer segments, consisting of the various 
regions of the given fragment of the basic library, do not contain some of the 
25 repeated k-1 sequences. In every segment a connecting sequence (a 
connecting subfragment) is generated for a certain pair of the subfragments 
from the given fragment. The process of ordering consists of three steps: (1) 
generation of the k- tuple contents of each segment; (2) generation of 
subfragments in each segment; and (3) connection of the subfragments of the 
30 segments. Primary segments are defined as significant intersections and 
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differences of k-tuple contents of a given fragment of the basic library with 
the k-tuple contents of the pools of the ordering library. Secondary (shorter) 
segments are defined as intersections and differences of the k-tuple contents 
of the primary segments. 
5 There is a problem of accumulating both folse positive and 

n^ative k^tuples in both the differences and intersections. The &lse n^ative 
k*tuples from starting sequences accumulate in the intersections (overlapping 
parts)» as well as &Ise positive k-tuples occurring randomly in both sequences, 
but not in the relevant overlaying region. On the other hand, the majority 

10 of false pointives from either of the starting sequences is not taken up into 
int^sections. This is an example of the reduction of experimental errors from 
individual fragments by using information from fragments overlapping with 
them. The false k- tuples accumulate in the differences for another reason. 
The set of false negatives from the original sequences are enlarged for false 

15 positives from intersections and the set of false positives for those k-tuplcs 
which are not included in the intersection by error, i.e. are false negative in 
the intersection. If the starting sequences contain 10% false negative data, the 
primary and secondary intersections will contain 19% and 28% false negative 
k- tuples, respectively. On the other hand, a mathematical expectation of 77 

20 false positives may be predicted if the basic fragment and the pools have 
lengths of 500 bp and 10,000 bp, respectively. However, there is a possibility 
of recovering most of the "lost" k-tuples and of eliminating most of the false 
positive k-tuples. 

First, one has to determine a basic content of the k-tuples for 

25 a given segment as the intersection of a given pair of the k-tuple contents. 
This is followed by including all k- tuples of the starting k-tuple contents in 
the interseaton, which contain at one end k-1 and at the other end k-+ 
sequences which occur at the ends of two k-tuples of the basic set. This is 
done before generation of the differences thus preventing the accumulation of 

30 false positives in that process. Following that, the same type of enlargement 
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of k-tuple set is ai^lied to differences with the distinction that tite borrowing 
is from the intmections. All borrowed k-tuples are eliminated from the 
intersection files as false positives. 

The intersection, i.e. a set of common k-tuples, is defined for 
each pair (a basic fragment) X (a pool of ordering library). If the number of 
k*tiiples in the set is significant it is enlarged with the false negatives 
according to the described rule. The primary diffidence set is obtained by 
subtracting from a given basic fragment the obtained intersection set. The 
folse negative k-tuples are app^ded to the difference set by borrowing from 
the intersection set according to the described rule and, at the same time, 
removed from the intersection set as false positive k-tuples. When the basic 
fragment is longer than the pooled fragments, this difference can represent the 
two separate segments which somewhat reduces its utility in further steps. 
The primary segments are all generated intersections and differences of pairs 
(a basic fragment) X (a pool of ordering library) containing the significant 
number of k-tuples. K-tuple sets of secondary segments are obtained by 
comparison of k-tuple sets of all possible pairs of primary segments. The two 
differences are defined from each pair which produces the intersection with 
the significant number of k-tuples. The majority of available information 
from overlapped fragments is recovered in this step so that there is little to be 
gained from the third round of forming intersections and differences. 

(2) Generation of the subfragmenis of the segments is 
performed identically as described for the fragments of the basic library. 

(3) The method of connection of subfragments consists of 
sequentially determining the correctly linked pairs of subfragments among the 
subfragments from a given basic library fragment which have some overlapped 
ends. In the case of 4 relevant subfragments, two of which contain the same 
beginning and two having the same end, there are 4 different pairs of 
subfragments that can be cranected. In general 2 are correct and 2 are 
wrong. To find correct ones, the presence of the connecting sequences of 
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each pair is tested in the subfragments generated from all primary and 
secondary segments for a given basic fragment. The length and the position 
of the connecting sequence are chosen to avoid interfmnce with sequences 
which occur by chance. They are k+2 or longer* and include at least one 
S element 2 beside overlying sequence in both subfragments of a givoi pair. 
The connection is permitted only if the two connecting sequences are found 
and the remaining two do not exist. The two linked subfragments replace 
former subfragments in the file and the process is cyclically repeated. 

Repeated sequences are generated in this step. This means that 

10 some subfragments are included in linked subfragments more than once. They 
will be recognized by finding the relevant connecting sequence which engages 
one subftagment in connection with two different subfragments. 

The recognition of misconnected subfragments generated in the 
processes of building pSFs and merging pSFs into longer subfragments is 

15 based on testing whether the sequencer of subfragments from a given basic 
fragment exist in the sequences of subfragments generated in the segments for 
the fragment. The sequences from an incorrectly connected position will not 
be found indicating the misconnected subfragments. 

Beside the described three steps in ordering of subfragments 

20 some additional control steps or steps applicable to specific sequences will be 
necessary for the generation of more complete sequence without mistakes. 

The determination of which subfragment belongs to which 
segment is performed b comparison of contents of k-iuples in segments and 
subfragments. Because of the errors in the k-tuple contents (due to the 

25 primary error in pools and statistical errors due to the frequency of 
occurrences of k-tuples) the exact partitioning of subfragments is impossible. 
Thus, instead of "all or none" partition, the chance of coming from the given 
segment (P(sf,s)) is determined for each subfragment. This possibility is the 
function of the lengths of k-tuples, the lengths of subfragments. the lengths of 
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fragments of oidering library, the size of the pool, and of the percentage of 
&lse k-tuples in the file: 

P(sf,s)=(Ck-F)/Lsf, 
where Lsf is the length of subfiragmenl, Ck is the number of common k-tuples 
5 ftir a given subfragment/segmoit pair, and F is the param^r that includes 
relations between lengths of k-tuples, firagments of basic library, the size of 
the pool, and the enor percentage. 

Subftagments attributed to a particular segment are treated as 
redundant short pSFs and are submitted to a process of unambiguous 
10 connection. The definition of unambiguous connection is slightly different in 
this case, since it is based on a probability that subfragmcnts with overiapping 
cnd(s) belong to the segment considered. Besides, die accuracy of 
unambiguous connection is controlled by following the connection of these 
subftagments in other segments. After the connection in different segments, 
15 all of the obtained subfragmcnts are merged together, shoner subfragmcnts 
included within longer ones are eliminated, and the remaining ones are 
submitted to the ordinary connecting process. If the sequence is not 
regenerated completely, the process of partition and connection of 
subftagments is repeated with the same or less severe criterions of probability 
20 of belonging to the particular segment, followed by unambiguous connection. 

Using severe criteria for defining unambiguous overlap, some 
information is not used. Instead of a complete sequence, several subfragmcnts 
that define a number of possibilities for a given fragment are obtained. Using 
less severe criteria an accurate and complete sequence is generated. In a 
25 certain number of situations, e.g. an erroneous connection, it is possible to 
generate a complete, but an incorrect sequence, or to generate "monster- 
subftagments with no connection among them. Thus, for each fragment of the 
basic library one obtains: a) several possible solutions where one is correct 
and b) the most probable correct solution. Also, in a very small number of 
30 cases, due to die mistake in the subfragmeni generation process or due to the 
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spedfic ratio of the probabilities of belongings no unambiguous solution is 
generated or one, the most probable solution. These cases remain as 
incomplete sequences, or the unambiguous solution is obtained by comparing 
these data with other, overlapped fragments of basic library. 
5 The described algorithm was tested on a randomly gen»ated, 

SO kb sequence, containing 40% GC to simulate the GC content of the human 
genome. In the middle part of this sequoice were inserted various All, and 
some other repetitive sequences, of a total length of about 4 kb. To simulate 
an in vitro SBH experiment, the following operations were performed to 

10 prepare ^ropriate data. 

• Positions of sixty 5 kb overlapping "clones'* were randomly 
defined, to simulate preparation of a basic library: 

- Positions of one thousand 500 bp "clones" were randomly 
determined to simulate making the ordering library. These fragments were 

IS extracted from the sequence. Random pools of 20 fragments were made, and 
k-tuple sets of pools were determined and stored on the hard disk. These data 
are used in the subfragment ordering phase: For the same density of clones 
4 million clones in basic library and 3 million clones in ordering library are 
used for the entire human genome. The total number of 7 million clones is 

20 several fold smaller than the number of clones a few kb long for random 
cloning of almost all of genomic DNA and sequencing by a gel-based method. 

From the data on the starts and ends of S kb fragments, 1 17 
"informative fragments" were determined to be in the sequence. This was 
followed by determination of sets of overlapping k-tuples of which the single 

25 "informative fragment" consist. Only the subset of k-tuples matching a 
predetermined list were used. The list contained 65% 8-mers, 30% 9-mers, 
and S% 10-12-mers. Processes of generation and the ordering of 
sub^gments were performed on these data. 

The testing of the algorithm was performed on the simulated 

30 data in two experiments. The sequence of SO informative fragments was 
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regenerated with the 100% conecl data set (over 20,000 bp), and 26 
infonnative firagments (about 10.000 bp) with 10% false k-toples (5% positive 
and 5% negative ones). 

In the first experiment, all subfiragnients were correct and in 
5 only one out of SO infonnative fragments the sequence was not completely 
legeneiated but remained in the form of 5 subiragments. The analysis of 
positions of overlapped fragments of ordering library has shown that they lack 
the information for Ote unique ordering of the 5 subiragments. The 
subftagm«its may be connected in two ways based on ov^pping ends, 
10 1-2-3-4 5 and M-3-2-5. The only differojce is the exchange of positions of 
subftagments 2 and 4. Since subfragments 2, 3, and 4 are relatively short 
(total of about 100 bp), the relatively greater chance existed, and occurred in 
this case, that none of the fragments of ordering library started or ended in the 
subfragment 3 region. 

To simulate real sequencing, some false ("hybridization") data * 
was included as inpul in a number of experiments. In oligomer hybridization 
experiments, under proposed conditions, the only situation producing 
- unreliable data is the end mismatch versus full match hybridization. 
Therefore, in simulation only those k-tuples differing in a single element on 
either end from the real one were considered to be false positives. These 
"false" sets are made as follows. On the original set of a k-tuples of the 
informative fragment, a subset of 5% false positive k-tuples are added. False 
positive k-tuples are made by randomly picking a k-iuple from the set, 
copying it and altering a nucleotide on its beginning or end. This is followed 
5 by subtraction of a subset of 5% randomly chosen k-tuples. In this way the 
statistically expected number of the most complicated cases is generated in 
which the correct k-mple is replaced with a k-tuple with the wrong base on the 
Old. 

Production of k-tuple sets as described leads to up to 10% of 
) false data. This value varies from case to case, due to the randomness of 
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choice of k-tuples to be copied, altered, and erased. Nevertheless, this 
pexcentage 3-4 tiroes exceeds the amount of unreliable data in real 
hybridization experiments. The introduced error of 10% leads to the two fold 
increase in the numb^ of subfiagments both in fragments of basic library 
S (basic library informative ftagments) and in segments. About 10% of the final 
subfiragments have a wrong base at the end as expected for the k-tuple set 
which contains false positives (see generation of primary subftagments). 
Neither the cases of misconnection of subfitagments nor subftagments with the 
wrong sequence were observed. In 4 inf<mnative fragments out of 26 

10 examined in the ordering process the complete sequence was not regenerated. 
In all 4 cases the sequence was obtained in the form of several longer 
subitagments and several shorter subfragments contained in the same segment. 
This result shows that the algorithmic principles allow working with a large 
percentage of false data. 

15 The success of the generation of the sequence from its k-tuple 

content may be described in terms of completeness and accuracy. In the 
process of generation, two particular situations can be defmed: 1) Some part 
of the information is missing in the generated sequence, but one knows where 
the ambiguities are and to which type they belong, and 2) the regenerated 

20 sequence that is obtained does not match the sequence from which the k- tuple 
content is generated, but the mistake can not be detected. Assuming the 
algorithm is developed to its theoretical limits, as in the use of the exact 
k-tuple sets, only the first situation can take place. There the incompleteness 
results in a certain number of subfragments that may not be ordered 

2S unambiguously and the problem of determination of the exact length of 
monotonous sequences, i.e. the number of perfect tandem repeats. 

With false k-tuples, incorrect sequence may be generated. The 
reason for mistakes does not lie in the shortcomings of the algorithm, but in 
the fact that a given content of k-tuples unambiguously represents the sequence 

30 that differs from the original one. One may deHne three classes of error, 
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depending on the kind of the fitlse k- tuples piesent in the file. False negative 
k-tuples (which are not accompanied with the iabe positives) produce 
"deletions". False positive k-tuples are produdng "elongations (unequal 
crossing over)". False positives accompanied with £dse negatives are the 
5 reason for generation of "insertions", alone or combined witii ^deletions'. 
The deletions are produced when all of Uie k*tuples (or dxeir m^ority) 
between two possible starts of tiie subfiragments are false n^atives. Since 
every portion in the sequence is defined by k k- tuples, the occurrence of the 
deletions in a common case requires k consecutive false negatives. (With 
10 10% of the false negatives and k=8, Uiis situation takes place after every 108 
elements); This situation is extremely mfrequent even in mammalian genome 
sequencing using random libraries containing ten genome equivalents. 

Elongation of the end of the sequence caused by false positive 
k-tuples is the special case of "insertions", since the end of the sequence can 
15 be considered as the endless linear array of false negative k-tuples. One may 
consider a group of false positive k-iuples producing subfragments longer tiian 
one k-tuple. Situations of tiiis kind may be detected if subfragments are 
generated in overlapped fragments, like random physical fragments of the 
ordering library. An insertion, or insertion in place of a deletion, can arise 
20 as a result of ^jecific combinations of false positive and false negative 
k-tuples. In the first case, the number of consecutive false negatives is 
smaller than k. Both cases require several overlapping false positive k-tuples. 
The insertions and deletions are mostly theoretical possibihties without sizable 
practical repercussions since Uie requirements in the number and specificity 
23 of false k-tuples are simply too high. 

In every other situation of no meeting the theoretical 
requirement of the minimal number an the kind of tiie false positive and/or 
n^atives, mistakes in die k-tuples content may produce only the lesser 
completeness of a generated sequence. 
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CLAIMS 

1 . A method for analyzing nucidc acids by hybridization, 
comprising the steps of: 

arraying a tint plurality of nudeic acid segments on a first 
5 sector of a substrate; 

disposing a second plurality of nucleic add segments on a 
second sector of said substrate; 

exposing, under conditions discriminating between full 
complementarity and a one baise mismatch, said first plurality of nucidc add 
10 segments to a first hybridization probe in said first sector, said first 
hybridization probe being shorter than one from among said first plurality of 
nucleic acid segments, to said plurality of nucleic acid segments; 

incubating under conditions discriminating between full 
complementarity and a one base mismatch, a second hybridization probe in 
IS said second sector, said second hybridization probe being shorter than a 
segment from among said second plurality of nucleic acid segments and said 
second hybridization probe bdng different in sequence from said first 
hybridization probe; 

detecting hybridization of a hybridization probe to a nucleic acid 

20 segment; and 

analyzing the result. 

2. The method as recited in claim 1, further comprising, 
prior to said disposing step, the step of introducing a barrier to movement of 
a nucleic acid. 

25 3. The method as recited in claim 1 further comprising, 

after said arraying and said disposing step but before said incubating step, the 
step of introducing a barrier to movement of a nucleic acid. 
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4. The method as recited in claim 3 wherein said 
introducing step comprises pres^g a physical barrier against said substrate. 

5. The method as recited in claim 2 wherein said 
introducing step conqmses the step of applying a direction-switching electrical 
field perpendicular to said support to prevent the mixing of probes between 
sectors. 

6. The method as recited in claim 3 wherein said 
introducing step comprises the step of flying a direction-switching electrical 
field perpendicular to said support to prevent the mixing of probes between 
sectors. 

7. The method as recited in claim 1 wherein said arraying 
step comprises the step of spotting nucleic acid samples by means of a jnn 
array. 

8. The method as recited in claim 1 wherein said arraying 
step comprises the step of dispensing nucleic acid samples by an array of 
tubes. 

9. The method as recited in claim I wherein said arraying 
step comprises the step of jet printing nucleic acid samples. 

10. The method as recited in claim 1 wherein said exposing 
step comprises the step of applying a plurality of contiguously hybridizing 
probes. 
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1 1 . The method as recited in claim 1 wherein said incubating 
step comprises the ssq) of applying a plurality of contiguously hybridizing 
probes. 

12. The method as recited in claim 10 further comprising the 
5 stq9 of ligating at least two of said plurality of contiguously hybridizing 

probes. . 

13. The method as recited in claim 1 1 further comprising the 
step of ligating at least two of said plurality of contiguously hybridizing 
probes. 

10 14. The method as recited in claim 1 wherein said exposing 

step comprises the step of applying a plurality of competitively hybridizing 
probes having overlapping nucleic acid sequences. 

15. The method as recited in claim 1 wherein said incubating 
step comprises the step of applying a plurality of competitively hybridizing 

IS probes having overlapping nucleic acid sequences. 

16. The method as recited in claim I wherein a least two of 
said first plurality of nucleic acid segments are arrayed as a mixture. 

17. The method as recited in claim 1 wherein a least two of 
said second plurality of nucleic acid segments are disposed as a mixture. 

20 18. The method as recited in claim 1 further comprising the 

steps of preparing samples by digestion with an Hga 1 type restriction enzyme 
and ligating the resulting restriction fragments with an anchor. 
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19. The method as recited in claim 1 further comprising the 
step of selecting probes from a univnsal set of probes of a given length. 

20. The method as recited in daim 1 further comprising the 
stq> of selecting probes from an incomplete set of probes of a given length. 

5 21. The method as recited in claim 1 further comprising die 

step of selecting deoxyribonucleotide probes. 

22. The mediod as recited in claim 1 further comprising the 
step of selecting ribonucleotide probes. 

23. The method as recited in claim 1 further comprising the 
10 step of selecting a nucleic add analog selected from the group consisting of 

protein nucleic add probes and probes containing base analogs. 

24. The method as rcdted in claim 1 further comprising the 
step of multiplex labelling of probes. 

25. The method as recited in claim 1 further comprising the ' 
15 step of degrading a label on an unhybridized probe. 

26. The method as recited in claim 19 wherein said exposing 
or said incubating step comprises the step of assembling a set of universal 
probes 6, 7, 8, 9 or 10 bases in length. 

27. The method as recited in claim 19 wherein said exposing 
20 or said incubating step comprises the step of assembling a set of universal 

probes 6, 7. 8, 9 or 10 bases in length. 
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28. The method as recited in claim 20 whmin said exposing 
or said incubating step comprises the stq> of assembling an incomplete sel of 
probes 5, 6. 7. 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 
24, 25, 26, 27, 28, 29 or 30 bases in length. 

5 29. Apparatus analyzing nucleic adds by hybridization 

comprising a substrate having points of attachment for nucleic acid fragments, 
said substrate being segmented by hydrophobic regions. 

30. The method as recited in daim 20 wherein said 
disposing step comprises the step of assembling an incomplete set of probes 
10 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23. 24, 25, 
26, 27, 28, 29 or 30 bases in length. 

3K The method of claim 1 further comprising the step of 
confirming the relative order of at least two bases in a segment by detecting 
hybridization of two or more probes having overlapping nucleic acid 
15 sequences including said at least two bases. 

32. A method for nucleotide sequence analysis comprising 

the steps of: 

introdudng a sample to an array of probes; 
adjusting the temperature to be one at which a majority of 
20 sample molecules are unassociated with ligated probes at any given time; 
adding a labelled probe to the mixture; 
incubating the mixture with ligase; 
removing free probes; and 
detecting ligation products. 
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33. The method as recited in claim 1 iiirth^ comprising the 
stq>s of defining additional probes for improving a desired result and 
repeating said exposing, incubating, detecting and analyzing steps. 

34. The method as recited in claim 1 further comprising the 
S step of stripping the substrate of probes for reuse of said pluralities of nucleic 

acid s^ments. 
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